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Although  a  number  of  mechanisms  have  been  proposed  to  explain 
the  existence  of  impurity  gradients  in  epitaxial  GaAs,  the  impurity 
gradient  model  presented  in  this  report  shows  that  such 
effects  are  inherent  to  the  growth  process.  That  is,  when  the 
carrier  concentrations  in  the  substrate,  growing  layer,  and  at  the 
growth  surface  are  different  at  the  growth  temperature,  non-uniform 
time-dependent  electric  fields  are  obtained  in  the  layer.  T^se— ^ 


fiel*ds  can  enhance  or  retard  the  motion  of  ionized  impurities  and 
defects  during  growth,  producing  impurity  gradients  at  the  outer 
surface  and  different  conductivity  regions  at  the  epitaxy- 
substrate  interface. 

High  temperature  resistivity  and  Hall  coefficient  measure¬ 
ments  were  made  on  epitaxial  layers  and  substrates  and  analyzed 
using  a  self-consistent  four-band  model  to  obtain  quantitative 
results  for  the  impurity  gradient  model  for  growth  on  heavily- 
doped  n-type  substrates  predicts  that  thin  (0.05  pm  or  larger) 
p-type  regions  at  the  epitaxy-substrate  interface  will  be 
produced  under  conditions  commonly  encountered  in  the  epitaxial 
growth  of  GaAs . 

t  s. 

Zero-bias  capacitance  measurements  on  Schottky  barrier  diode 
reveal  the  presence  of  heavily-doped  n-type  regions  near  the 
outer  surface.  These  surface  gradients  caused  by  surface  states 
are  predicted  by  the  model.  The  effects  of  surface  states  can 
also  explain  the  observed  self-compensation  in  GaAs. 
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1 .  INTRODUCTION 


Since  epitaxial  layers  of  GaAs  were  first  used  [1]  for 
the  fabrication  of  Gunn  oscillators,  the  occurrence  of 
different  conductivity  regions  at  the  layer-substrate 
interface  [2,3]  has  been  a  recurring  problem  in  the  fab¬ 
rication  of  GciAs  microwave  devices.  To  explain  these 
interfacial  regions  a  number  of  mechanisms  have  been  pro¬ 
posed  including  outdiffusion  or  autodoping  of  amphoteric 
impurities  from  the  substrate  [2] ,  defects  formed  by  tran¬ 
sient  effects  during  the  initial  stages  of  growth  [4] , 
copper  contamination  of  the  substrate  [5] ,  substrate  de¬ 
composition  prior  to  growth  [6,7],  and  a  large  concentra¬ 
tion  of  silicon  at  the  interface  from  the  quartz  reactor  [8] . 
The  impurity  gradient  model  presented  in  this  report  shows 
that  such  interfacial  regions  are  inherent  to  the  growth 
process.  That  is,  if  all  other  possible  mechanisms  for 
the  formation  of  interfacial  regions  were  eliminated,  in 
most  epitaxial  growth  situations  of  practical  interest 
these  regions  would  still  be  observed. 

The  types  of  interfacial  regions  observed  for  growth 
on  heavily-doped  substrates  are  either  high-resistivity 
n-type  or  p-type  regions  with  thicknesses  anywhere  from 
about  5  ym  to  less  than  0.1  ym  [2,4-13].  It  has  been 

*  The  numbers  in  parentheses  in  the  text  indicate 
references  in  the  Bibliography. 


shown,  however,  that  the  thicker  high-resistance  n-type 
regions  are  more  properly  interpreted  as  thin  p-type 
regions  [12,13].  The  impurity  gradient  model  presented 
in. this  report  predicts  a  thin  (0-05  ym  or  larger)  p-region 
at  the  layer-svibstrate  interface  ur  -‘tx  conditions  commonly 
encountered  in  epitaxial  growth.  Figure  1.1,  from  the  data 
of  Sato  and  lido  [10],  shows  the  type  of  interfacial  region 
predicted  by  the  model.  At  the  interface  between  the 
n-type  epitaxial  layer  and  the  n"*"  substrate,  there  is  a 
p-region  approximately  0.05  ym  thick.  The  observed  profile 
is  due  to  the  outdiffusion  of  electrons  and  the  fact  that 
difterential  capacitance  measurements  measure  the  carrier 
concentration,  not  the  net  impurity  concentration  [14] . 

The  types  of  interfacial  regions  commonly  observed  for 
growth  on  chromium-doped  substrates  are  both  p-type  [15]  and 
n-type  [3].  Figure  1.2,  from  the  data  of  Yamasaki  [16], 
shows  the  type  of  heavier  doped  n-region  predicted  by  the 
impurity  gradient  model.  The  profile  in  Figure  1.2  is  due 
to  the  outdiffusion  of  electrons  from  a  0.05  ym  heavier-doped 
region.  The  thickness  of  this  inter facial  region  is  very 
close  to  that  predicted  by  the  impurity  gradient  model. 

This  model  can  also  explain  the  screening,  above  the 
intrinsic  carrier  concentration,  of  the  donor-acceptor 
impurity  interaction  during  incorporation  into  the  epitax¬ 
ial  layer.  That  is,  when  the  impurity  concentrations  are 
sufficiently  below  the  intrinsic  concentration,  the 


DISTANCE  FROM  SURFACE  (/*m) 


Figure  1.1  Carrier  concentration  versus  distance  from  the 
surface  for  an  n-type  epitaxial  layer  grovm  on 
an  n^  substrate  with  a  p-type  interfacial  region 


Epitaxial  Layer  Thickness  (pm) 


Figure  1.2  Carrier  concentration  versus  epitaxial  layer 
thickness  for  an  n-type  epitaxial  layer 
grown  m  a  Cr-doped  substrate  with  an 
n^  interfacial  region. 
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impurities  are  screened  by  the  intrinsic  concentration. 
This  produces  a  constant  compensation  ratio,  N^/N^. 
However,  when  the  impurity  concentrations  are  above  the 
intrinsic  concentration,  the  donors  and  acceptors  are 
expected  to  interact  causing  the  compensation  to  increase. 
The  dashed  curve  in  Figure  1.3  is  for  the  equilibrium  con¬ 
dition  where  there  is  no  variation  between  the  electron 
concentration  at  the  surface  and  inside  the  epitaxial 
layer.  As  can  be  seen,  this  dashed  curve  (no  surface 
states)  does  not  agree  with  the  experimental  data.  On 
the  other  hand,  band  bending  [17,18]  at  the  surface,  due 
to  a  lower  electron  concentration  at  the  surface  caused 
by  surface  states,  can  screen  the  donor-acceptor  impurity 
interaction  at  higher  concentrations.  The  solid  curve  in 

Figure  1.3  was  fit  to  the  experimental  data  by  assuming 

11  -2 

the  presence  of  6.3  x  10  cm  surface  states. 

The  impurity  gradient  model  can  also  explain  the 
chromium  pile-up  observed  at  the  surface  of  annealed, 
chromium-doped  substrates  [19,20,21]  and  epitaxial  layers 
grown  on  chromium-doped  substrates  [22,23].  Since  chro¬ 
mium  has  a  charge  of  +2  in  its  compensating  state  and  +3 
in  its  neutral  state  [24],  the  electric  field  in  the  sur¬ 
face  region  will  cause  the  chromium  atoms  to  drift 
towards  the  surface.  This  drift  will  produce  a  chromium 

gradient  at  the  surface  with  a  thickness  on  the  order 
of  0.10  ym. 


Nj*Ng  (cm’’) 


Figure  1 . 


3  Net  donor  concentration  versus  total 

ionized  impurity  concentration  (N^+N^) •  The 

dashed  line  is  calculated  for  zero  surface 

states.  The  solid  line  is  calculated  for 
1 1  —2 

6.3  X  10  cm  surface  states.  The  layers 
are  doped  with  the  column- IV  impurities: 

Si  (•)  ,  Ge  (A)  ,  and  Sn  (■)  . 
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The  salient  features  of  the  impurity  gradient  model 
for  two  sets  of  growth  conditions  are  shown  in  Figure  1.4. 

Figure  1.4  (a)  is  for  the  growth  of  a  lightly  doped 
epitaxial  layer,  which  is  intrinsic  at  1000  K,  on  a 
heavily-doped  substrate.  At  1000  K  the  substrate  is  still 
extrinsic  and  the  resulting  difference  in  electron  concen¬ 
tration  produces  an  electric  field  between  the  substrate 
and  the  epitaxial  layer.  There  is  a  similar  electric 
field  between  the  epitaxial  layer  and  the  surface  caused 
by  the  lower  electron  concentration  at  the  surface  due 
to  surface  states.  These  electric  fields  will  cause  the 
positively-charged  donors  to  drift  away  from  the  substrate- 
epitaxial  layer  interface  and  towards  the  surface. 
Similarly,  the  negatively-charged  acceptors  will  drift 
away  from  the  surface  and  towards  the  substrate-layer 
interface.  This  will  produce  a  donor  pile-up  at  the  sur¬ 
face  and  an  acceptor  pile-up  at  the  interface. 

If  this  profile  is  frozen-in  when  the  growth  process 
is  terminated,  at  300  K  a  p-region  will  be  formed  at  the 
substrate- layer  interface.  At  the  surface  there  will 
still  be  band  bending  due  to  surface  states  but  the  band 
bending  will  be  reduced  due  to  the  heavier-doped  region 
produced  by  the  donor  impurity  gradient. 

Figure  1.4  (b)  indicates  an  epitaxial  layer,  which 
is  extrinsic  at  1000  K,  being  grown  on  a  chromium-doped 


Figure  1.4  Energy  band  diagreuns  for  (a)  a  heavily-doped 
substrate  and  (b)  a  Cr-doped  substrate  at  the 
growth  temperature  (1000  K)  and  300K.  The 
impurity  gradients  formed  at  1000  K  are  assumed 
to  be  frozen  in  when  the  layer  is  cooled. 
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siibstrate,  which  is  either  slightly  p-type  or  intrinsic  at 
the  growth  temperature.  The  difference  in  electron  con¬ 
centration  will  produce  a  net  electric  field  between  the 
epitaxial  layer  and  the  substrate.  There  will  also  be 
band  bending  at  the  surface  due  to  surface  states  similar 
to  that  previously  discussed.  This  band  bending  will  re¬ 
sult  in  an  electric  field  which  will  cause  positively- 
charged  donors  to  drift  towards  the  surface  and 
negatively-charged  acceptors  to  drift  away.  In  a  similar 
manner,  the  electric  field  at  the  interface  will  cause  the 
donors  to  drift  towards  the  interface  and  the  acceptors 
away. 

If  this  profile  is  frozen  in  when  the  growth  process 
is  terminated,  at  300  K  heavier-doped  regions  will  result 
at  both  the  surface  and  siabstrate-epitaxial  layer  inter¬ 
face.  Thus,  the  band  bending  at  the  surface  will  be 
reduced  due  to  these  impurity  gradients. 

This  impurity  gradient  model  predicts  a  surface 

impurity  gradient  which  has  not  previously  been  reported 

in  the  literature.  This  surface  gradient  should  occur 

in  all  n-type  epitaxial  layers  which  are  doped  below 
18  •S 

approximately  10  cm 

In  this  report,  high  temperature  Hall  and  resistivity 
measurements  on  epitaxial  layers  and  substrates  will  be 
discussed  in  Section  2.  From  these  data  the  carrier 
concentrations  of  the  epitaxial  layers  and  sxibstrates 
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can  be  calculated.  Using  these  results  a  quantitative 
impurity  gradient  model  will  be  discussed  in  Section  3. 
From  the  model  a  surface  impurity  gradient  is  predicted. 
In  Section  4.1  the  surface  state  density  needed  to  form 
the  experimentally  observed  surface  impurity  gradients 
will  be  calculated.  This  value  will  then  be  compared 
to  the  surface  state  density  necessary  to  produce  the 
experimentally  measured  compensation  seen  in  GaAs.  As 
will  be  seen,  the  good  agreement  between  these  two 
independently-obtained  values  tends  to  confirm  the 
impurity  gradient  model. 


2 .  HIGH-TEMPERATURE  CARRIER  TRANSPORT 


2.1  EPITAXIAL  LAYERS 

High  temperatxire  Hall  measurements  are  important, 
not  only  for  the  insights  they  give  into  band  structure 
parcuneters,  but  also  because  of  the  data  they  provide 
for  modeling  epitaxial  growth  phenomena  [25]  and  other 
high  temperature  processing.  Earlier  high  temperature 
Hall  measurements  were  made  on  polycrystalline  [26]  or 
bulk  grown  GeiAs  [27-30]  .  To  improve  the  quality  of  these 
data,  more  recent  measurements  have  been  made  on  higher- 
purity  epitaxial  GeiAs  layers  [31-r35]  .  For  room- temperature 
electrical  measurements,  epitaxial  layers  are  grown  on 
Cr-doped  substrates  with  a  room  temperature  resistivity 

Q 

of  approximately  10  n-cm.  Since  the  purest  n-type  layers 
have  a  resistivity  of  only  about  50  fi-cm,  the  conduction 
through  the  substrate  is  negligible.  This  substrate,  how¬ 
ever,  provides  severe  restrictions  to  high- temperature 
electrical  measurements  on  epitaxial  material  due  to  a 
rapid  increase  in  the  conductivity  of  the  substrate. 

Zucca  [36]  has  measured  a  decrease  in  resistivity  of 
lO^n-cm  when  Cr-doped  substrates  were  heated  from  300  to 
500  K.  Allen  [37]  has  observed  an  increase  of  10^  in  the 
carrier  concentration  of  Cr-doped  GaAs  when  heated  from 
300  to  600  K,  thereby  reducing  the  resistivity  to  only 


i 


-12- 


2 

10  n-cm.  This  substrate  conduction  can  make  the  parallel 
current  paths  of  the  epitaxial  layer  and  substrate  of 
comparable  resistance,  thus  strongly  affecting  the  measured 
electrical  properties  of  epitaxial  layers  above  500  K. 

Much  of  the  recent  data  on  epitaxial  layers  have  been 
obtained  with  the  substrates  still  attached  [31  33] . 
Although  anomalous  features  observed  have  been  attributed 
to  deep  donors,  some  of  these  results  are  probably  due 
to  conduction  in  the  substrate. 

Akita  et  at.  [34]  have  briefly  reported  measurements 
on  one  epitaxial  sample  from  which  the  substrate  was  re¬ 
moved,  while  Blood  [35]  has  made  measurements  on  two  more 
epitaxial  layers.  The  data  of  Akita  et  al.  [34]  and 
Blood  [35]  have  been  limited  to  the  temperature  range 
300  to  800  K  due  to  evaporation  of  arsenic.  The  only 
electrical  measurements  made  above  900  K  on  GaAs  have 
been  obtained  by  Roberts  [30]  and  Smith  [38].  Roberts' 
and  Smith's  measurements  were  made  with  a  maximum  tempera¬ 
ture  between  1000  and  1190°C  on  bulk  grown  crystals  in 
am  arsenic  atmosphere.  In  addition,  all  of  these  measure¬ 
ments  have  been  analyzed  using  a  F-X-L  ordering  of  the 
conduction  band  minima. 

In  this  section,  electrical  measurements  are  re¬ 
ported  in  the  range  300-1250  K  on  epitaxial  GaAs  from 
which  the  substrate  has  been  removed  by  mechanical-chemical 
polishing.  The  epitaxial  layers  have  room  temperature 


-13- 


carrier  concentrations  ramging  from  10^^  to  10 ^^cm”^. 

In  Section  2.1.1  the  preparation  of  these  layers  and  de¬ 
tails  of  the  measurements  are  described,  while  in 
Section  2.1.2  the  data  are  analyzed  using  the  nev  r-L-X 
ordering  of  the  conduction  band  minima. 

2.1.1  Experiment 

2.1. 1.1  Scunple  Preparation 

The  epitaxial  layers  were  grown  on  Cr-doped  semi- 

insulating  GaAs  substrates  in  an  AsCl2“Ga-H2  vapor  phase 

reactor  [39,40].  The  substrates  had  a  crystallographic 

orientation  of  <211>  Ga  to  achieve,  simultaneously,  a  high 

growth  rate  [41]  and  high  purity  [42] .  With  a  background 

14  -3 

concentration  around  10  cm  ,  the  doping  of  the  epitaxial 
layers  was  controlled  by  the  addition  of  tin  to  the  gallium 
melt.  In  this  manner,  the  epitaxial  layers  were  grown  to 
a  thickness  of  between  280  and  360  pm  with  doping  concentra¬ 
tions  in  the  lO^^-lO^^cm”^  range.  The  as-grown  crystals 

were  cleaved  into  diamond-shaped  Scunples  with  areas  of 
2 

about  1  cm  .  The  samples  were  then  mounted  face  down  on 
a  lapping  block  and  five  pm  alumina  was  used  to  remove 
the  substrate  to  within  approximately  100  pm  of  the 
epitaxial  layer.  The  remaining  substrate  was  then  re¬ 
moved  with  a  mechanical-chemical  polish  using  1%  bromine- 
methanol.  This  polishing  was  continued  10  to  20  pm 
into  the  epitaxial  layer  to  eliminate  any  interfacial 
regions . 
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To  prevent  the  evaporation  of  arsenic  at  high  tem¬ 
peratures,  the  epitaxial  layers  were  encapsulated.  The 
Seunples  were  first  etched  in  H2S0^ :H202 :H20 (5 :1 ;1)  to 

remove  approximately  10  microns,  then  Si^N^  was  deposited 

0 

in  a  pyrolytic  reactor  143]  to  a  thickness  of  3000  A. 

This  Si3N^  encapsulation  worked  well  below  1150  K.  Above 
this  temperature,  however,  the  Si^N^  cracked  causing  a 
loss  of  arsenic.  To  prevent  this,  a  composite  layer  of 

arsenic-doped  Si02  and  Si^N^  first  described  by  Lindow 

0 

et  al.  [44]  was  used.  First  1000  A  of  pyrolytic  Si^N^  was 
deposited;  then  approximately  1  pm  of  arsenic-doped  Si02 
was  pyrolytically  deposited  on  top  of  the  Si^N^.  This 
layer  maintained  its  integrity  up  to  1254  K,  which  was  the 
maximum  temperature  used  in  our  measurements.  Mask  in¬ 
tegrity  was  confirmed  when  Hall  and  resistivity  measure¬ 
ments  with  increasing  temperature  were  the  Seune  (within 
experimental  error)  as  with  decreasing  temperature  and 

subsequent  high  temperature  measurements.  Contact  windows 

.  2 
with  a  total  area  of  about  1/16  cm  were  provided  by 

masking  growth  of  the  encapsulating  material  on  the  four 

corners  of  the  upper  surface. 

2. 1.1. 2  Measurement  Apparatus 

The  apparatus  shown  in  Figure  2 . 1  was  used  for  the 

high- temperature  resistivity  and  Hall  measurements.  To 

regulate  the  temperature  this  apparatus  was  placed  inside 

an  open-tube,  two-zone  resistance  furnace.  The  furnace 


Carbon 


tube  and  the  silica  portions  of  the  apparatus  were  made  of 
high-purity  quartz.  The  windings  of  the  furnace  were  made 
of  platinum  with  an  outer  casing  of  aluminum  to  minimize 
distortion  of  the  magnetic  field.  The  magnetic  field  was 
provided  by  placing  the  furnace  between  the  poles  of  a 
4-inch  electromagnet.  The  magnetic  field  for  all  measure¬ 
ments  was  set  at  3.5  kG  and  was  measured  with  a  gaussmeter 
to  be  within  5%  of  this  value  over  the  4  inches  of  the 
furnace  inside  the  poles  of  the  magnet.  Hydrogen,  purified 
with  a  palladium  diffusion  cell,  was  used  as  the  ambient 
gas  to  prevent  oxidation  of  the  sample  or  probing  apparatus. 
The  temperature  was  measured  with  a  platinum-10 %  rhodium 
thermocouple  which,  with  a  digital  thermometer,  provided 
an  accuracy  of  ±  0.4°C.  The  probing  apparatus  was  centered 
in  a  two-inch  region  of  constant  (±  0.4°C)  temperature 
which  had  been  centered  in  the  magnetic  field. 

Above  650  K  the  carbon  probes  shown  in  Figure  2.1 
were  used  to  provide  ohmic  contacts  for  the  resistivity 
and  Hall  measurements.  The  carbon  probes  were  pressure 
contacts  with  a  sharp  edge  for  contacting  the  GeiAs  surface. 
These  probes  were  held  rigid  by  carbon  nuts  and  screws 
mounted  through  a  quartz  base  plate  which  isolated  the 
four  contacts.  This  plate  was  mounted  on  a  four-bore 
quartz  capillary  tvibe  to  provide  insulated  electrical 
paths  outside  the  high-temperature  zone.  High  purity 
molybdenum  wire  was  inserted  in  the  capillaries  of  the 
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tubing  and  connected  to  the  carbon  probes  by  a  pressure 
contact  provided  by  another  carbon  nut  on  the  carbon  screw. 
The  15-mil  dieuneter  molybdenum  wire  was  connected  outside 
the  furnace  to  a  switching  network  arranged  to  perform 
van  der  Pauw  measurements.  To  prevent  diffusion  of  gases 
through  the  capillaries  of  the  tubing,  the  low  temperature 
end  was  sealed. 

Below  650  K  the  carbon  probes  give  rectifying  rather 
than  ohmic  contacts .  To  provide  an  ohmic  contact  at  these 
temperatures,  a  AuSn  eutectic  alloy  was  electroplated  in 
the  windows  of  the  encapsulant  and  alloyed  at  450°C  in  a 
reducing  atmosphere.  This  produced  ohmic  contacts  with  the 
carbon  probes  from  300  to  650  K. 

2. 1.1. 3  Experimental  Results 

Several  properties  of  the  measured  samples  are  summar¬ 
ized  in  Table  2.1.  The  300  K  carrier  concentrations  were 
determined  from  the  measured  Hall  coefficient  assuming  a 
Hall  factor  of  unity.  The  total  impurity  concentrations, 
Np+N^,  were  calculated  from  the  77  K  mobility  [45]  for 
Samples  Nimiber  2  and  3.  Since  this  technique  is  not  valid 
for  heavier-doped  samples,  similar  data  were  obtained  from 
the  300  K  mobility  [46]  for  Sample  Number  1. 

The  resistivity,  p,  was  measured  by  the  van  der  Pauw  [47] 
technique.  All  the  samples  had  a  resistance  ratio  between 
1.05  and  1.00  in  both  the  low  (300-600  K)  and  high  (650- 
1250  K)  temperature  ranges.  This  ratio  was  also  temperature 


Table  2.1  Sample  Properties 


Sample 

No. 

"3OOK 

(cm“^) 

NT^+N, 

D  A 

(cm“^) 

■1 

Thickness 

(vm) 

1 

1.2x10^’ 

2.0x10^"^ 

0.25 

340 

2 

9.7x10^^ 

1.6x10^® 

0.23 

310 

3 

15 

1.9x10-^^ 

15 

4.2x10-^^ 

0.38 

260 

independent.  The  resistances  were  calculated  from  the  slope 
of  the  current- voltage  characteristics  by  linear  regression. 
These  characteristics  were  obtained  from  voltage  measure¬ 
ments  for  at  least  three  different  currents  for  each 
measurement.  The  voltage  was  measured  using  a  high  impe¬ 
dance  digital  multimeter  to  four  digit  accuracy,  while 
the  current  was  supplied  by  a  constant  current  source. 

The  measured  I-V  characteristics  were  linear:  that  is, 
the  coefficient  for  determination  of  linear  regression 
was  0.9999 ±0.0001,  while  an  optimum  fit  would  be  1.0.  The 
resistivity  showed  the  expected  increase  from  the  room 
temperature  value  with  rising  temperature  due  to  increased 
scattering  of  the  carriers.  As  intrinsic  conduction  began 
to  dominate,  the  resistivity  started  decreasing  with  in¬ 
creasing  temperature  due  to  the  increase  in  the  carrier 
concentration . 

The  Hall  coefficient,  Rjj,  was  measured  in  a  magnetic 
field  of  3.5  kG  with  the  van  der  Pauw  [47]  technique. 

Since  the  seimples  were  diamond-shaped,  the  current  was 
directed  through  the  longest  path  to  give  the  best  ac¬ 
curacy.  We  estimate  the  cumulative  errors  of  the  measure¬ 
ments  under  these  conditions  to  be  under  5%.  The  single¬ 
band  carrier  concentration  was  calculated  from  the  Hall 
coefficient  by. 


where  q  is  the  electronic  charge,  and  the  Hall  scattering 
factor  was  assumed  to  be  unity.  This  assumption  introduces 
a  small  temperature-dependent  error  in  the  analysis  [48]  . 
The  temperatvire  dependencies  of  this  single-band  carrier 
concentration  for  the  seimples  indicated  in  Table  2 . 1  are 
plotted  in  Figure  2.2.  Since  all  the  donors  are  thermally 
ionized  at  room  temperature,  the  carrier  concentration  re¬ 
mains  constant  until  sufficient  carriers  are  thermally  ex¬ 
cited  into  higher  conduction  band  minima.  The  slight 
apparent  rise  in  the  carrier  concentration  up  to  400  K 
is  due  to  the  decrease  in  the  Hall  factor  for  polar  mode 
scattering  [48] .  As  the  samples  are  heated  into  the  600  K 
range,  there  is  a  noticable  rise  in  the  Hall  coefficient 
due  to  the  thermal  excitation  of  carriers  into  the  lower 
mobility  Lg  minima.  The  sharp  rise  of  the  single-band 
carrier  concentration  above  750  K  for  the  highest-doped 
sample  is  due  to  the  onset  of  intrinsic  conduction.  This 
is  confirmed  by  the  higher  temperatures  needed  to  produce 
equivalent  increases  for  the  heavier-doped  samples. 

The  measured  mobility,  Uj^^,  is  determined  as  the  Hall 
coefficient  divided  by  the  resistivity:  ^/P-  These 
mobility  values  for  the  temperature  range  300-1250  K  are 
plotted  in  Figure  2.3.  The  mobility  variations  of  the 
lighter-doped  samples  have  characteristics  similar  to  the 
Scunples  measured  by  Blood  [35]  . 


TEMPERATURE  (K) 


Figure  2.2  Reciprocal  of  Hall  coefficient  times  the 

elementary  charge  as  a  function  of  temperature 
for  three  GaAs  epitaxial  samples. 


TEMPERATURE  (K) 


Figure  2.3  Measured  mobility,  Rjj/P/  versus  temperature  for 
three  GaAs  epitaxial  seunples.  The  solid  line 
is  from  the  analysis  by  Rode  [49,50]. 
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In  the  range  300-500  K  the  mobility  can  be  expressed 
as  with  N  =  1.25,  while  in  the  range  500-700  K 

the  temperature  dependence  is  slower.  Above  750  K  the 
measured  mobility  values  for  the  lighter-doped  seunples 
again  fall  rapidly  with  temperature.  These  variations 
with  temperature  are  similar  to  those  calculated  by 
Rode  [49,50]  for  pure  GeU^s.  In  comparison,  however,  our 
S2unples  have  a  slightly  reduced  slope  between  300  and  500  K 
due  to  ionized  impurity  scattering  [51].  Thus,  the 
9.7  X  10^^  cm  ^  doped  Sample  Number  2  has  a  smaller  slope 
than  the  1.9  x  10^^  cm~^  doped  Seunple  Number  3. 

Above  500  K  the  slope  of  the  data  for  the  lighter- 
doped  samples  decreases.  This  feature  is  predicted  by 
Rode's  analysis  [49,50].  However,  the  calculation  does 
not  account  for  the  increase  in  slope  above  750  K.  This 
behavior  can  be  partially  ascribed  to  thermal  transfer 
of  electrons  from  the  high  mobility  Tg  minimum  to  the  low 
mobility  Lg  minima  [50].  The  heaviest  doped  Sample  Number  1, 
however,  does  not  show  this  decrease  at  750  K,  but  rather 
at  950  K.  This  is  partially  due  to  the  high  ionized- 
impurity  scattering  of  Sample  Number  1,  combined  with  the 
fact  that,  when  intrinsic  conduction  dominates,  a  different 
scattering  mechanism  begins  to  cause  a  rapid  decrease  of 
the  mobility  in  the  Fg  minimum.  This  mechanism  is  the 
scattering  of  carriers  from  the  Fg  minimum  to  the  Lg 
minima,  then  back  to  the  Fg  minimum  [52].  It  is  obvious 
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that  this  mechanism  is  strongly  dependent  on  the  number  of 
carriers  in  the  upper  portion  of  the  Tg  minimum.  There¬ 
fore,  both  temperature  and  carrier  concentration  are  im¬ 
portant  factors.  The  combined  effect  of  this  scattering 
mechanism  and  the  thermal  transfer  to  lower-mobility  upper 
minima  determines  the  mobility  at  high  temperature.  Thus, 
above  1100  K  the  measured  mobility  data  for  all  of  the 
samples  begin  to  approach  the  same  value. 

To  summarize  the  findings  of  this  section,  we  find 
good  agreement  among  our  experimental  temperature  depen¬ 
dence  of  the  measured  mobility.  Blood's  experimental  data, 
and  Rode's  calculations. 

2.1.2  Analysis 

2. 1.2.1  Energy  Separations 

In  contrast  to  previous  conduction-band  ordering 
in  GaAs,  Aspnes  et  at,  [53]  have  recently  found,  from 
Schottky  barrier  electroreflectance  measurements ,  that  the 
Lg  conduction  band  minima  are  170±30  meV  below  the  Xg  mini¬ 
ma.  Aspnes  [54,55]  has  also  shown  that  all  present  experi¬ 
mental  data  can  be  reinterpreted  using  Tg-Lg-Xg  ordering  of 
the  conduction  band  minima.  The  temperature  dependence  of 
these  minima  with  respect  to  the  valence  band  is  needed  to 
determine  the  intrinsic  carrier  concentration  and  relative 
occupancy  of  various  minima  at  high  temperatures.  The  top 
of  the  valence  band,  Fg,  will  be  used  as  the  reference 
energy  {E=0)  throughout  the  following  discussion. 


-25- 


The  Tg-Tg  separation  has  been  measured  with  great 
accuracy  by  Sell  et  al.  [56,57]  and  by  Panish  and  Casey  [58], 
These  data  were  used  by  Thurmond  [59]  to  determine  the 
coefficients  E^,  a,  and  6  of  the  equation  given  by 
Varshni  [60], 

Ep  =  Ep  -aTV(T  +  8)  ,  (2.1) 

0 

as  1.519  eV,  5.405  x  10  ^eV/K,  and  204  K,  respectively. 

These  coefficients  predict  energies  that  agree  with  ex¬ 
periment  with  a  mean-square  deviation  of  2.6  meV. 

The  Xg-Fg  separation  can  be  determined  from  the  above 
data  and  data  on  indirect  transitions  from  the  electrons, 
at  Fg  to  higher  conduction  band  minima  in  degenerate  n-type 
material.  Using  the  data  of  Onton  [61] ,  Aspnes  determined 
the  Xg-Fg  separation  by  lineshape  analysis  to  be  0.462  ± 

0.005  eV  at  2  K.  This  results  in  a  separation  of  1.981  eV 
at  2  K.  The  coefficient  8  in  Equation  (2.1)  is  related  to 
the  Debye  temperature  for  the  fundamental  absorption  edges 
of  Si,  Ge,  GaP  and  GaAs  [54].  From  data  for  Si  and  GaP, 
Aspnes  assximes  that  8  =  204  K  for  the  Xg-Fg  separation  in 
GeU^s.  By  analyzing  transport  data  and  combining  the  re¬ 
sults  with  a  rigid-band  hypothesis  [62]  and  results  of 
measurements  on  GaP  [63],  Aspnes  [54]  calculates  the  a  co¬ 
efficient  of  Equation  (2.1)  for  Xg  to  be  4.6  x  lO”^  eV/K. 

The  Xg-Lg  energy  separation  has  been  determined  from 
core-level  electroreflectance  measurements  to  be  0.17  ± 

0.03  eV  at  110  K  [53].  To  determine  the  coefficient  a  for 
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Teible  2.2  Energies  and  temperature  dependence  of  the 

c  c  ^ 

Lf  and  X,  conduction-band  minima  of  GaAs  relative 

to  the  top  of  the  valence  band,  Fg.  The 
coefficients  E^,  a  and  3  are  defined  by  Equa¬ 
tion  2.1. 


Eq  (ev) 

a(eV/K) 

6  (K) 

1.519 

5.405x10"^ 

204 

1.818 

6.5x10"^ 

204 

1.981 

4.6x10"'* 

204 
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the  I'g“rg  separation  Aspnes  [54]  uses  the  Auvergne  et  al .  [62] 
postulate  of  rigid  valence  bands  and  his  measured  energy 
shift  of  115  meV  of  the  L^-Lg  transition  from  4  to  295  K 
from  Schottky  barrier  electroreflectance  measurements  in 
the  3  eV  range  [64].  From  this  he  calculates  a  =  6.5  x 
10  ^  eV/K  for  the  Lg"rg  separation  [54]  .  Using  this  a  and 
the  x|-Lg  energy  separation  at  110  K,  the  Lg-Tg  separation 
at  0  K  is  calculated  to  be  1.818  eV.  All  of  the  coefficients 
for  the  temperature  dependencies  of  the  valence  band-to- 
conduction  band  energies  are  summarized  in  Table  2.2. 

2. 1.2. 2  Four-Band  Model 

To  calculate  the  true  carrier  concentration  a  four- 
band  model  is  needed.  The  Hall  coefficient  under  these 
conditions  is 


^  '3<“r“r*"L  L^x^x-PPh'^  ' 

where  n^,  n^^  and  n^  are  the  electron  concentrations;  Up, 

Ut  and  y  are  the  electron  mobilities;  p  is  the  hole  con- 
centration;  and  y  is  the  hole  mobility  in  the  valence  band 
maxima.  The  corresponding  expression  for  the  resistivity  is 


1 


P 


(2.3) 


To  eliminate  two  unknowns  in  these  equations  the  relative 

occupancy  of  the  conduction  band  minima  can  be  calculated 

for  nondegenerate  material  using  Boltzman  statistics. 

Although  Sample  Number  1,  with  a  doping  level  of  1.2  x 
17  -3  . 

10  cm  IS  degenerate  at  room  temperature,  beyond  450  K 
(the  temperature  range  of  interest)  it  is  non-degenerate. 
Defining  nj^/np=a  and  n^/np=b,  the  relative  occupancy  can 
be  expressed  as. 


(2.4) 


where  Ep,  and  E^  are  the  respective  separations  between 
the  various  minima  and  Pg.  Wp »  and  m^^  are  the  respective 
effective  masses  of  the  minima. 

2.1.2. 3  Effective  Masses 

For  a  parabolic  band  the  density-of-states  mass  can 
be  expressed  as 


(2.5) 


where  m  and  m,  are  the  transverse  and  longitudinal 
masses,  and  is  the  number  of  equivalent  minima.  The 
isotropic  Fg  mass,  mp,  at  2  K  is  m^  =  m^  =0. 0665m  [65] 

The  variation  of  this  mass  with  temperature  can  be 
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calculated  from  k*p  theory  [66]  according  to 


m  1  +  Ep 

where  Ep  is  an  energy  related  to  the  momentum  matrix  ele¬ 
ment,  and  A  =  0.341  eV  is  the  spin-orbit  splitting, 
rg"r7  [64],  From  the  temperature  dependence  of  the  T 
energy  gap  and  the  above  data,  Ep  is  calculated  to  be 
7.51  eV.  Thus,  using  Equation  (2.6)  the  temperature  de¬ 
pendence  of  the  effective  mass  in  the  rf  minimum  is 


(2.6) 


m 


1  +  7.51 


where  the  temperature  dependence  of  E 
tion  (2.1) , 


1 

r 

r 


- ^  .  (2.7) 

■+irr34ij 

is  given  by  Equa- 


Aspnes  [54]  calculates  the  transverse  effective  mass 
for  the  Lg  minima  from  his  Schottky  barrier  electro¬ 
reflectance  measurements  [64]  to  be  0.0754m.  He  estimates 
the  longitudinal  mass  at  Lg  by  comparison  of  GaAs  with  Ge. 
His  estimated  longitudinal  mass  value  of  1.9m  leads  to  a 
density-of-states  effective  mass  of  m^  =  0.56m  at  0  K  [54]. 
The  temperature  dependence  of  this  density-of-states  effec¬ 
tive  mass  is  determined  by  the  temperature  variation  of 

C  V 

the  Lg  -  Lg  interband  energy  separation  which  affects  the 

transverse  mass,  m.  [54].  Using  Equation  (2.5),  the  data 

^L 


of  Table  2.1,  k*p  theory  and  the  values  calculated  by 
Aspnes,  the  density-of-states  effective  mass  for  the  L® 
minima  as  a  function  at  temperature  can  be  expressed  as 


=  {5.51  [1  +  19.3 
m 


1 

(ev) 


(ev) 


+  0.22 


)]> 


2/3 


(2.8) 


where  is  the  bg-Lg  interband  energy  separation. 

We  use  the  Xg  density-of-states  effective  mass  of 
0.85m  calculated  by  Pitt  and  Lees  [67]  from  their  high- 
pressure  Hall  effect  and  resistivity  data.  No  temperature 
dependence  is  assumed  for  the  effective  mass  at  Xg.  This 
is  a  reasonable  assumption  since  even  at  1250  K  we  esti¬ 
mate  the  percentage  of  carriers  in  the  Xg  minima  to  be  at 
most  only  5%. 

2. 1.2. 4  Hall  Coefficient 

2. 1.2. 4.1  Extrinsic  Conductivity 

If  the  material  has  a  carrier  concentration  well 
above  the  intrinsic  concentration,  then  the  contribution 
of  the  holes  can  be  neglected.  Below  1000  K  the  carriers 
in  the  Xg  minima  represent  less  than  3%  of  the  total  number 
of  carriers  in  the  conduction  bands.  Since  the  mobility 
of  the  carriers  in  the  Xg  minima  is  the  lowest  of  the  con¬ 
duction  band  minima  at  F,  L,  and  X,  the  total  effect  of 
the  carriers  in  Xg  can  be  neglected  for  the  Hall  coefficient 
[Equation  (2.2)]  and  resistivity  equations  [Equation  (2.3)]. 


» 
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Neglecting  the  carriers  in  the  Tg  maxima  and  Xg  minima, 
the  expression  for  the  Hall  coefficient  from  Equation  (2.2) 
becomes 


q  Rh  =  —  2 

nj,{l+ag)^ 


(2.9) 


where  a  is  defined  by  Equation  (2.4)  and  B  =  Pj^/Pp. 

Using  the  temperature  dependence  of  the  energy  gaps 
given  in  Table  2.2  and  the  effective  masses  given  in 
Equations  (2.7)  and  (2.8),  a  can  be  calculated.  Combining 
this  value  of  a  and  the  measured  Hall  coefficient  into 
Equation  (2.9)  yields  an  equation  that  has  two  unknowns; 
np  and  6.  Another  equation  can  be  obtained  if  the  total 
concentration  of  the  layer  is  assumed  to  be  the  doping 
concentration,  Nj^-N^.  This  is  a  valid  assumption  since 
there  are  no  appreciable  deep  donor  levels  and  the  intrinsic 
concentration  is  small  enough  to  neglect.  Thus,  it  is  as¬ 
sumed  that 


Up  +  Ul  -  n^QQ 


(2.10) 


The  400  K  carrier  concentration  was  chosen  because  the 
Hall  factor  for  polar  mode  scattering  is  a  minimiam  at  this 
temperature 


Combining  Equations  (2.9)  and  (2.10)  yields  the 


equation, 


q  Rjj  = 


(1+a)  (1+aB^) 
''400K  (l+a0)2 


(2.11) 


which  has  only  one  unknown:  g.  This  mobility  ratio,  g, 
was  calculated  for  the  heavier-doped  Samples  Nximber  2 
and  3.  The  results  for  these  epitaxial  layers  are  plotted 
in  Figure  2.4.  The  measvirements  were  done  until  intrinsic 
conduction  could  no  longer  be  neglected  in  Equation  (2.9). 
The  final  point  calculated  for  Saunple  Number  1  used  a 
modified  Equation  (2.10).  That  is,  the  intrinsic  carrier 
concentration  was  added  to  the  400  K  concentration  in 
Equation  (2.10).  However,  because  of  the  small  mobility 
of  holes,  the  intrinsic  conduction  can  still  be  neglected 
in  Equation  (2.9).  The  temperature  dependence  shown  in 

Figure  2.4  indicates  that  the  assumption  of  a  constant 
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mobility  ratio  for  samples  doped  above  10  cm  is  a 
poor  assumption.  It  will  be  shown  later  in  this  paper 
that  the  results  in  Figure  2.4  also  apply  to  Sample  Num¬ 
ber  3  doped  at  1.9  x  10^^  cm 

For  pure  GaAs  intravalley  polar-optical  scattering 
predominates  in  the  Fg  conduction  band  minimum  [50]  until 
intrinsic  conduction  becomes  predominant.  Even  the  purest 
samples  measured  by  Blood  [35]  and  our  Sample  Number  3 
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Figure  2 


. 4  Ratio  of  the  mobility  in  the  Lg  minima  to  the 
mobility  in  the  Fg  minimum  as  a  function  of 
temperature  for  Samples  No.  1  (O)  and  No.  2  (a) . 
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show  a  smaller  temperature  dependence  [49]  between  300- 
500  K  than  would  be  expected  for  pure  GaAs.  This  is  due 
to  ionized- impurity  scattering,  which  will  be  important 
for  the  mobility  of  the  carriers  in  r®.  However,  it  will 
have  a  much  smaller  effect  on  the  carriers  in  Lg.  The 
dominant  scattering  mechanism  in  the  Xg  minima  is  equiva¬ 
lent  intervalley  scattering  [68,69]  and  this  mechanism 
can  likewise  be  assiimed  to  dominate  the  mobility  in  the 
Lg  minima.  Although  scattering  between  equivalent  valleys 
leads  to  a  larger  temperature  coefficient  than  measured 
in  our  samples  [70]  ,  in  GaP  the  observed  value  [7i]  is  not 
significantly  different  from  that  of  Tg.  If  the  same 
applies  to  GclAs,  then  for  pure  GaAs  the  mobility  ratio, 

B,  will  remain  constant  with  increasing  temperature.  How¬ 
ever,  ionized  impurity  scattering  reduces  the  temperature 
dependence  of  the  Fg  minimum,  causing  the  observed  tempera¬ 
ture  dependence  of  the  mobility  ratio,  6,  given  in  Figure  4 
2. 1.2. 4. 2  Intrinsic  Conductivity 

Near  the  point  where  intrinsic  conduction  begins  to 

become  important,  the  ratio  of  the  hole  mobility  to  the 

c 

electron  mobility  in  Fg  can  be  calculated.  This  was  done 
by  using  the  intrinsic  mass  action  law, 

(np  +  nj^)p  =  n^^,  (2.1 

obtained  by  neglecting  the  contributions  of  the  carriers  in 
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the  Xg  minima.  The  previous  Equations  (2.9)  and  (2.10) 
become , 


”400K  ”i' 


Using  the  energy  gap  parameters  in  Table  2.2  and  the  effec¬ 
tive  masses  given  in  Equations  (2.7)  and  (2.8),  the  intrin¬ 
sic  carrier  concentration  was  calculated.  This  intrinsic 
concentration  was  used  in  Equations  (2.12)  and  (2.13)  and 
the  mobility  ratio  Uj^/Up  was  calculated  in  the  900-1000  K 
range  for  Sample  Number  1  and  in  the  750-850  K  range  for 
Sample  Number  2.  It  was  found  that  this  ratio  was  approx¬ 
imately  0.053,  which  is  the  same  as  the  room  temperature 
value.  Thus,  the  ratio  Mj^/Up  was  assumed  to  be  constant. 

2. 1.2. 5  Resistivity 

Neglecting  the  contribution  of  the  holes  and  the  elec¬ 
trons  in  x|.  Equation  (2.3)  becomes. 


P 


_ 1 _ 

qUp(nj,  +  nj^g) 


(2.14) 


With  this  equation,  the  8  calculated  in  Figure  2.4,  the 
measured  resistivity,  and  Equations  (2.4)  and  (2.10),  the 


-36- 


mobillties  and  were  calculated.  These  mobilities 
are  plotted  versus  temperature  in  Figure  2.5,  along  with 
the  hole  mobility. 

With  these  mobilities  a  three-band  carrier  concentra¬ 
tion  can  be  calculated.  To  calculate  the  desired  four- 
band  carrier  concentration  an  assumption  must  be  made 

about  the  temperature  dependence  of  the  Xg  mobility.  We 

c  c 

have  assumed  that  the  Xg  mobility  is  one-third  the  Lg 

mobility,  based  on  the  higher  density  of  states  at  Xg  and 
high  pressure  transport  measurements  [54] .  We  have  also 
assumed  that  the  Xg  and  Lg  mobilities  have  the  same  tempera¬ 
ture  dependence ,  since  intravalley  scattering  dominates 
the  mobility  for  both  sets  of  minima.  Using  these  as¬ 
sumptions  and  the  mobilities  in  Figure  2.5,  the  four-band 
carrier  concentration  was  calculated  from  Equations  (2.1), 
(2.2),  (2.4),  (2.7),  (2.8),  and  the  data  in  Table  2.2. 

When  intrinsic  conduction  begins  to  dominate  as  previously 

c  c 

discussed,  intervalley  scattering  between  Tg  and  Lg  be¬ 
gins  to  control  the  Tg  mobility.  This  causes  the  mobility 
ratio,  g,  to  flatten  out  or  to  start  increasing  [52].  Un¬ 
der  extrinsic  conditions  the  best  results  for  the  four-band 
carrier  concentration  were  obtained  when  the  mobility  ratio 
g ,  was  taken  to  have  the  value  in  Figure  2.4.  When  Scunple 
Numbers  1,  2,  and  3  reached  1050,  950,  and  900  K,  respect¬ 
ively,  the  mobility  ratio  leveled  off  and  began  to  rise 
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Figure  2.5  Mobilities  in  theP®  minimvun,  Lgininima  and  Pg 

maxima  as  a  function  of  temperature  for  Samples 
No.  1  (O)  and  No.  2  (Q)  . 
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slightly.  This  rise  was  approximately  5%  for  every  100  K 
increase  in  temperature.  The  final  four-band  carrier  con¬ 
centration  versus  temperature  is  plotted  in  Figure  2.6. 

The  calculated  intrinsic  carrier  concentration  was  calcu-; 
lated  using  the  data  in  Table  2.2  and  Equations  (2.1), 

(2.7)  and  (2.8). 

2.1.3  Conclusions 

From  the  preceding  results  and  the  other  data  dis¬ 
cussed,  a  self-consistent  model  has  been  developed  for  the 
evaluation  of  carrier  transport  in  n-type  GaAs  at  elevated 
temperatures.  This  model  supports  the  r|-Lg-Xg  ordering 
of  the  conduction-band  minima  demonstrated  by  Aspnes  and 
co-workers.  The  temperature  dependencies  of  the  various 
energy  separations  have  been  shown,  within  experimental 
error ,  to  have  the  variations  given  in  Table  2.2. 

Using  these  temperature  dependencies  of  the  energy 
separations  and  the  effective  mass  variations,  a  four- 
band  intrinsic  carrier  concentration  has  been  calculated 
which  is  consistent  with  measured  data  at  high  tempera¬ 
tures.  We,  therefore,  expect  this  calculated  intrinsic 
carrier  concentration  to  be  correct  within  experimental 
error. 

The  temperature  dependence  of  the  mobility  in  the 
Lg  minima  was  shown  to  decrease  faster  with  increasing 
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Figure  2 


.6  Four-band  carrier  concentration  as  a  function  of 
temperature  for  GeiAs  epiteocial  Seunples  No.  1  (o)  r 
No.  2  (□)  and  No.  3  (A) .  The  dashed  line  is  the 
calculated  intrinsic  carrier  concentration. 


temperature  than  the  mobility  in  the  Tg  minimum.  This 
different  temperature  dependence  is  attributed  to  differ- 
ent  scattering  mechanisms.  In  the  Lg  minima  equivalent 
intervalley  scattering  dominates  the  mobility,  while  in 
Tg  a  combination  of  polar-mode  and  ionized-impurity  intra¬ 
valley  scattering  dominates.  Therefore,  taking  the  ratio 

of  the  mobility  in  Lg  to  the  mobility  in  Pg  as  constant  is 

15  -3 

not  applicable  to  layers  doped  above  10  cm  .  When  in¬ 
trinsic  conduction  dominates,  intervalley  scattering  be¬ 
tween  Tg  ouid  Lg  will  cause  the  mobility  in  Pg  to  have  a 
slightly  greater  temperature  dependence  than  the  mobility 
in  L°. 

Combining  the  temperature  dependence  of  the  mobilities 
and  energy  gaps  with  the  four-band  Hall  coefficient  equa¬ 
tion  produced  results  which  were  consistent  with  the  calcu¬ 
lated  intrinsic  carrier  concentration  and  the  measured 
carrier  concentration. 

2 . 2  SUBSTRATES 

To  model  epitaxial  growth  phenomena  the  carrier  con¬ 
centrations  of  the  substrates  at  high  temperature  must 
be  known.  Thus,  high- temperature  Hall  measurements  were 
performed  on  both  heavily-doped  and  chromium-doped 
substrates.  The  substrates  were  prepared  using  the  en¬ 
capsulation  procedure  of  Section  2. 1.1.1.  The  measure¬ 
ment  apparatus  was  exactly  the  same  as  described  in 
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Section  2. 1.1. 2,  however,  no  measurements  were  made 
below  650  K. 

The  resistivity,  p,  was  measured  by  the  van  der  Pauw 
technique  [47] .  All  the  samples  had  a  resistance  ratio 
between  1.05  and  1.00  which  was  temperature  independent. 
The  resistances  were  calculated  from  the  slope  of  the 
current-voltage  characteristics  by  linear  regression. 

These  characteristics  were  obtained  from  voltage  measure¬ 
ments  for  at  least  three  different  currents  for  each 
measurement.  The  voltage  was  measured  using  a  high 
impedance  digital  multimeter  to  four  digit  accuracy, 
while  the  current  was  supplied  by  a  constant  current 
source.  The  measured  I-V  characteristics  were  linear; 
that  is,  the  coefficient  for  determination  of  linear 
regression  was  0.9999  ±  0.0001,  while  an  optimum  fit  would 
be  1.0. 


The  Hall  coefficient,  R^j,  was  measured  in  a  magnetic 
field  of  3.5  kG  with  the  van  der  Pauw  technique  [47]. 

Since  some  of  the  samples  were  diamond- shaped,  the  current 
was  directed  through  the  longest  path  to  give  the  best 
accuracy.  The  single-band  carrier  concentration  was 
calculated  from  the  Hall  coefficient  by. 


where  q  is  the  electronic  charge. 


and  the  Hall  scattering  fac¬ 


tor  was  assumed  to  be  unity.  This  assumption  introduces  a  small 
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temperature-dependent  error  in  the  analysis  [48] .  The 
temperature  dependencies  of  this  single-band  carrier 
concentration  for  the  chromium-doped  substrates  are 
plotted  in  Figure.  2.7. 

The  measured  mobility,  is  determined  as  the  Hall 

coefficient  divided  by  the  resistivity:  Rg/P*  These 

mobility  values  for  the  temperature  range  650-1225  K  are 

plotted  in  Figure  2.8.  The  range  of  mobility  values  from 

2  2 

approximately  2000  cm  /V-sec  at  650  K  to  700  cm  /V-sec  at 
1200  K  is  the  same  range  exhibited  by  the  heaviest  doped 
seunple  in  Figure  2.3.  This  similarity  indicates  that 
there  is  a  large  aunount  of  ionized  impurity  scattering  in 
the  chromium-doped  substrates. 

The  data  in  Figure  2.7  and  the  Hall  measurements  of 
a  heavily-doped  substrate  at  1000  K  were  interpreted  using 
the  four-band  model  developed  in  Section  2.1.  The  multi¬ 
band  electron  concentrations  for  the  chromium- doped 
substrates  as  a  function  of  temperature  are  shown  in 
Figure  2.9.  The  dashed  line  is  the  calculated  intrinsic 
carrier  concentration.  From  these  data,  below  1100  K  the 
chromium-doped  substrates  are  slightly  p-type,  while  above 
1100  K  they  appear  to  be  intrinsic.  This  indicates  that 
these  substrates  do  not  have  the  pinned  Fermi  level  ob¬ 
served  by  Zucca  [36]. 


TEMPERATURE  (K) 


Figure  2.8  Measured  mobility,  Rjj/p,  versus  temperature  for 
three  Cr-doped  GaAs  substrates. 
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Figure 


2 . 9  Four-band  electron  concentration  as  a  function 

of  temperature  for  three  Cr-doped  GaAs  substrates. 
The  dashed  line  is  the  calculated  intrinsic 
carrier  concentration. 
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The  heavily-doped  siibstrate  was  only  measured  at 

room  temperature  and  1000  K.  Using  the  multi-band  model, 

18  —3 

the  sxibstrate  was  still  heavily-doped  (1.8  x  10  cm  )  at 
1000  K.  Thus,  as  expected  a  heavily-doped  substrate  will 
remain  extrinsic  until  temperatures  far  exceeding  those 
used  for  epitaxial  growth  are  reached. 

In  conclusion,  chromium- doped  substrates  are  slight¬ 
ly  p-type  while  heavily-doped  substrates  are  n^  at  the 
temperatures  commonly  encoxantered  for  vapor-phase 
epitaxial  growth. 


3.  IMPURITY  GRADIENT  MODEL 

From  the  high-temperature  Hall  data  of  Section  2, 
it  can  be  seen  that  the  carrier  concentrations,  at  typical 
growth  temperatures,  in  the  substrate  and  the  growing 
epitaxial  layer  are  different  for  many  commonly-encountered 
growth  situations.  The  situations  presented  in  Figure  1.4, 
for  the  growth  of  a  lightly-doped  epitaxial  layer  on  a 
heavily-doped  substrate  and  the  growth  of  an  extrinsic 
layer  on  a  nearly  intrinsic  chromium-doped  substrate,  are 
two  examples  where  this  occurs.  Once  carrier  concentra¬ 
tions  in  the  substrates  and  the  growing  layers  are  known 
from  the  high- temperature  Hall  data,  quantitative  impurity 
gradient  models  can  be  developed  for  these  and  other 
epitaxial  growth  situations. 

3.1  TIME-DEPENDENT  ANALYSIS 

3.1.1  Electric  Field 

For  a  time-dependent  analysis  of  the  problem  only 
two  assumptions  are  required  to  find  a  formal  solution 
for  the  electric  field.  These  assumptions  are  as  follows; 

(1)  The  epitaxial  growth  rates  are  small  enough 
that  the  physical  processes  are  approximately  in  thermo¬ 
dynamic  equilibri\im. 

(2)  Since  the  donor  and  acceptor  mobilities,  and 
u  ,  are  much  less  than  the  electron  and  hole  mobilities, 

3 
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Ujj  and  Pp,  it  is  assumed  that  the  ionized  donor  and  ac¬ 
ceptor  concentrations,  and  N^,  are  sufficiently  small 

that  Va  "^n'  P^'p* 

With  these  assiunptions  the  equations  to  be  solved 
are  Poisson's  equation. 


3E 


=  _2_  z  N. 
eeo  j  3 


(3.1) 


vrtiere  N.  =  +N, ,  -N^,  +p,  -n ; 

j  G  a 

the  conduction  current  equations  for  ionized  donors, 
acceptors,  electrons,  and  holes, 

J  =  qy.  N.E  r  qD.  ^  ?  (3.2) 

J  J  J  ^  J  3x 

and  the  continuity  equations. 


3N  .  3J  . 

^  ^ 


(3.3) 


where  the  upper  -  sign  in  the  current  .and  continuity  equa¬ 
tions  is  for  the  positively  charged  donors  and  holes, 
while  the  lower  +  sign  is  for  the  negatively  charged 
acceptors  and  electrons.  At  low  growth  rates,  the  total 
current  is  zero. 


J  =  z;j.  +  ee  H=o; 
j  3  0  3t  ' 


(3.4) 


and  for  dilute  solutions  the  mass  action  laws  are 


np  =  n.  ",  nN^  =  ,  pN^  =  K^N. 


(3.5) 


J 
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where  N,°  and  N  °  are  the  neutral  donors  and  acceptors, 
d  a 

and  n.^,  K,  and  are  the  equilibrium  constants.  In 

X  Cl  A 

these  equations  N ^ ,  J ^  and  E  are  assumed  to  be  functions 
of  time,  t,  and  distance  from  the  substrate,  -x. 

Putting  Equation  (3.2)  into  Equation  (3.4)  produces 


qE 

+ 


^a**a  ^  ^p^  ^n"^ 


q{-D,  +  D.  rr: - + 


d  3x  "a  3x 


P  3X 


3E 

^^0  TE  “  °  ' 


D  2l^) 
n  ax' 


(3.6) 


which  with  Assumption  2  and  Einstein's  relation  reduces  to 

P  »  3Sl  (S  “Sx  ~  ‘^n  3x)  ^^0 

q  (ypP  +  U„nj  ■  q(UpP  +  y„n) 

Using  the  mass  action  law  for  the  generation  of  holes  and 
electrons  given  in  Equation  (3.5),  Equation  (3.7)  becomes 


E 


q 


3n 


nEE„ 

_ 0 


3E 

3t  • 


Equation  (3.8)  can  be  rewritten  as 


(3.8) 


3E 

3t 


+ 


IcT  .Iv  3Q 
q  'n'  3x  7t  ' 


where 


+  Uj^n^)dt'. 


(3.9) 
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The  differential  equation.  Equation  (3.9) ,  can  be 
solved  by  using  the  integrating  factor  e®  where  Q  is 
given  by  Equation  (3.10).  Formally,  the  solution  for 
the  electric  field  from  Equation  (3.9)  is 


E(x,t) 


e  ^E(x,o)  = 


-IcT 

q 


-Q  t 
’  / 


t  Q 
e 


3Q 

3t 


_1^ 

n 


dt', (3.11) 


where  Q  is  given  by  Equation  (3.10). 

For  heavily-doped  siibstrates  where  the  surface  states 
are  compensated  by  the  stobstrate  doping,  E(x,o)  is 
approximately  zero  for  all  x.  With  this  condition 
Equation  (3.11)  can  be  written  as 


E 


V. 

/ 


It 


n 


(3.12) 


0 

For  substrates  which  are  near  intrinsic  at  the  growth  tem¬ 
perature,  such  as  chromium-doped  substrates,  E(x,o)  would 
not  be  zero  and  the  full  form  of  Equation  (3.11)  must  be 
used. 


3.1.2  Impurity  Concentrations 

Combining  Equations  (3.2)  and  (3.3)  for  electrons 

gives 


q 


In  -  i_  r 

3t  ~  3x  [ 


qy^nE  +  qD 


3n  ■ 

n  3x  J  , 


(3.13) 
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which  can  be  written  as 
1  3n  3^n  ,  ct  3 

5“  Tt  “  rr  w  ^ 

n  o  X 


(nE)  . 


(3.14) 


Using  the  value  for  the  electric  field  given  in 
Equation  (3.11)  in  Equation  (3.14),  the  equation  for 

determining  the  electron  concentration  can  be  obtained. 
This  equation  is 


3n  ^  3^  3  f 

°n  ^  ■  3x2  ‘  ^  t 


Q  3Q 
Tt^ 


n  3^ 


(3.15) 
+  i  E(x,ol]} 


where  Q  is  given  by  Equation  (3.10) . 

For  heavily-doped  substrates.  Equation  (3.15)  becomes 


_1  3n 

D'  3t 
n 


“  I?’ 


1  In 
n  3x 


(3.16) 


The  donor  and  acceptor  concentrations  can  then  be 
determined  from 


D  . 
3 


+ 


1 

n 


(3.17) 


3n  , 


qE  ( X ,  0 ) 
kT 


]} 


where  n  is  determined  from  Equation  (3.16),  Q  is  given 

by  Equation  (3.10),  and  the  upper  +  sign  is  for  the 
positively-charged  donors,  while  the  lower  -  sign  is 
for  the  negatively-charged  acceptors. 
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Although  in  general  the  time-dependent  solution  of 
Equations  (3.10),  (3.11),  (3.15),  and  (3.17)  must  be 
used,  at  low  growth  rates,  v,  and  thin  epitaxial  layer 
thicknesses,  L,  the  ionized  donors  and  acceptors  are  in 
equilibrium  with  the  electric  field  and  the  time- 
independent  solution  is  a  good  approximation.  That  is, 
the  time- independent  solution  can  be  used  for  v  L  much 
less  than  D ^ ,  the  impurity  diffusion  coefficients. 

3.2  TIME- INDEPENDENT  SOLUTIONS 

3.2.1  Heavily-Doped  Substrate 

3. 2. 1.1  Assumptions 

For  the  time- independent  impurity  gradient  model 
four  assumptions  are  used  which  provide  a  reasonably 
good  approximation  to  the  growth  of  an  n-layer  on  a 
heavily-doped  substrate  and  which  are  sufficiently  con¬ 
strained  to  produce  analytical  result.  The  assumptions 
are  as  follows : 

(1)  The  carrier  concentration  of  the  substrate  is 
fixed  by  the  substrate  doping,  N^,  which  is  greater  than 
the  carrier  concentration  at  the  growth  temperature,  n^. 

(2)  The  carrier  concentration  of  the  growing  layer, 
n^^,  is  determined  by  the  growth  temperature  and  not  by 
the  ionized  donors,  N^,  and  acceptors,  N^,  incorporated 
during  the  growth  process.  That  is,  n^^  is  greater  than 


N  ,  and  N  . 

u  d 
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(3)  The  carrier  concentration  at  the  growth  surface, 
n^,  is  fixed  by  surface  states  and  is  less  than  n^^. 

(4)  The  product  of  growth  rate,  v,  and  epitaxial 
layer  thickness,  L,  is  much  less  than  the  impurity  diffu¬ 
sion  coefficients,  D ^ . 

Practically,  these  assumptions  approximate  the  growth  of 
a  very  thin,  lightly-doped  epitaxial  layer  on  a  heavily- 
doped  substrate  with  a  very  low  growth  rate. 

Since  the  surface  of  the  epitaxial  layer  is  moving 
due  to  the  growth  process,  a  moving  coordinate  system  is 
used.  The  transformation  to  the  moving  coordinate  system 
is  given  by 

C=x+vt=x+L,  (3.18) 

where  C  is  the  distance  from  the  growth  surface,  -  x  is 
the  distance  from  the  substrate,  v  is  the  growth  rate,  and 
L  is  the  layer  thickness. 

3. 2. 1.2  Electron  Concentration 
3. 2. 1.2.1  Solution 

Since  the  donor  and  acceptor  concentrations  are 
sufficiently  low,  the  carrier  concentration  of  the  growing 
layer  is  determined  by  the  growth  temperature,  not  by 
the  ionized  donors  and  acceptors.  From  Section  2. 1.2. 5 
the  epitaxial  layer  will  be  intrinsic  at  the  growth  tem¬ 
perature  for  a  net  donor  concentration,  less  than 


10^®  cm”^.  This  analysis  should  then  apply  for  layers 
doped  below  this  level. 

The  time-independent  equation  for  electrons  is 


3n 

Tt 


0 


(nE)  + 


(3.19) 


from  the  combination  of  the  conduction  equation. 
Equation  (3.2),  and  the  continuity  equation.  Equation 
(3.3).  Poisson's  equation.  Equation  (3.1)  for  the 
lightly  doped  case  is 


af  =  ^  (P  -  n) 


(3.20) 


Putting  Equation  (3.20)  into  Equation  (3.19)  yields, 


nq 


HR 


ee. 


(p  -  n)  + 


d^n 


n 


=  0 


d5' 


(3.21) 


For  the  time- independent  case  the  electron  conduction 
current,  at  most  would  be  a  constant.  Iii  this  analy¬ 
sis  for  equilibrium  the  electron  conduction  current  is 
zero.  Equation  (3.2)  for  electrons  can  then  be  written  as 


‘'n  ,dn. 

"IT 


(3.22) 


Putting  Equation  (3.22)  into  Equation  (3.21)  and  using 
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the  mass  action  law  for  the  generation  of  holes  and  elec¬ 
trons,  Equation  (3.5)  produces 


'^^n  ,  2  2.  _*°n  .dn.  ^ 

—  (ni  -  n  )  -  -5-  (g^) 


This  equation  can  be  written, 
in  the  form 


2 

+  ^  =  0  (3.23) 

using  Einstein's  relation. 


d^n  1  ,dn.^  -  o 


where 


1/2 


A  = 


and 


2_2 


1/2 


.CSH) 

We  kT  / 
0 


Equation  (3.24)  has  solutions 


n  =  I  tanh^  (  -|^  C  + 

and 

n  =  I  coth^  (  5  +  C2) . 


(3.24) 


C3.25) 
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These  solutions  can  be  rewritten  as 


n 


n.  tanh^  +  C, ) 
1  ^i  ^ 


(3.26) 


And 


n  »  n^  coth^  +  C2) 


(3.27) 


where 


1/2 


q^n. 


(3.28) 


3 . 2 . 1 . 2 . 2  Boundary  Conditions 

At  the  surface  of  the  growing  layer  the  electron 
concentration,  n^,  is  assumed  to  be  fixed  by  surface 
states  below  the  intrinsic  electron  concentration.  The 
solution  for  the  electron  concentration  near  the  surface 
is  given  by  Equation  (3.26).  At  the  surface  (5  =  0), 
the  constant,  can  be  calculated,  using  the  boundary 

condition  n(5  =  0)  =  n^,  as 

-1  % 

a  =  tanh  ^  (jp)  ,  (3.29) 

“i 

where  equals  a.  The  electron  concentration  near  the 
surface  can  be  written  from  Equation  (3.26)  as 


n(5) 


2 

n^  tanh 


(^  +  a). 


(3.30) 
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where  a  is  given  by  Equation  (3.29)  and  is  given  by 
Equation  (3.28) . 

Near  the  substrate-layer  interface  (5  =  L)  the 
electron  concentration  in  the  epitaxial  layer,  n^^,  is 


nj^(5)  =  n^  coth^  (^j^~  ^  +  C^) , 


(3.31) 


from  Equation  (3.27).  This  equation  can  be  rewritten  as 

(3.32) 


n^(x)  =  n^  coth^(^  +  C2)  , 


where  -x  is  the  distance  in  the  epitaxial  layer  from  the 
interface.  The  electron  concentration  in  the  substrate, 
n2,  can  be  approximated  using  Equation  (3.26)  as 


n2(x)  =  Ng  tanh^(^  +  C3) 

e 


(3.33) 


where  x  is  the  distance  from  the  interface  in  the  sub¬ 
strate  and 

n.  1/2 

-  L.  (^)  . 

s 

since  the  electron  concentration  is  continuous  at 
the  interface  (C  =  L) ,  the  relationship  between  C2  and 
can  be  found  using  Equations  (3.32)  and  (3.33)  as 


N  1/2 

coth  C2  =  (jj-;-)  tanh  C^. 

The  exact  values  of  C2  and  can  be  calculated  from 


(3.34) 
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the  single-valueness  of  the  electric  field  at  the  in¬ 
terface  (x  =  0).  This  will  be  done  in  Section  3. 2. 1.2.1. 
Using  the  values  calculated  in  that  section,  the  electron 
concentration  is  plotted  in  Figure  3.1  for  various 
substrate-intrinsic  concentration  ratios,  N  /n. .  In 

S  X 

this  figure  -x  is  the  distance  from  the  interface  into 
the  epitaxial  layer  and  x  is  the  distance  from  the  inter¬ 
face  into  the  substrate. 

3.2.1. 3  Electric  Field 
3. 2. 1.3.1  Interface 

The  equilibrium  electric  field  can  be  calculated 
from  Equation  (3.22).  Using  Einstein's  relation  this 
equation  can  be  written  as 


E(x) 


kT  ,1.  dn 
q  ^n^  dx  * 


(3.35) 


The  electric  field  in  the  epitaxial  layer  (E^^)  near  the 
interface  (x  =  0)  is  obtained  from  Equation  (3.35)  as 


Ei(x)  =  -  csch  2  (J^  +  C2)  , 

where  -x  is  the  distance  from  the  interface  into  the 
epitaxial  layer  and  is  given  by  Equation  (3.28). 

The  electric  field  in  the  substrate  (E2)  is  cal¬ 
culated  from  Poisson’s  equation. 


(3.36) 


dE- 
_ 2 

dx 


ss  (N3-n2), 
0 


(3.37) 


Figure  3.1  Electron  concentration  versus  distance  for 

various  substrate-intrinsic  concentration  ratios 
Ng/n^^.  The  interface  is  at  x=0  and  -x  is  the 
distance  into  the  epitaxial  layer  while  +x  is 
the  distance  into  the  substrate. 


where  n2  is  given  by  Equation  (3.33).  Inserting  Equa¬ 
tion  (3.33)  into  Equation  (3.37)  gives 


^^^2  ^^s  2  X 

dir  =  ^  C3)  ' 

0  e 


(3.38) 


n.  1/2 

where  =  (j^-)  L^. 

s 

By  integrating  Equation  (3.38)  and  taking  the  electric 
field  far  into  the  substrate  to  be  zero,  the  electric  field 
near  the  interface  (x  =  0)  is  obtained.  This  electric 
field  is 


IF  { 1  's’}  - 


(3.39) 


where  x  is  the  distance  from  the  interface  into  the  sub¬ 
strate. 

Since  the  electric  field  must  be  single-valued  at 
the  interface  (x  =  0) ,  the  electric  fields  in  the  sub¬ 
strate,  E2,  and  in  the  epitaxial  layer,  E,  are  equal  at 
that  point.  Combining  Equations  (3.36)  and  (3.39)  and 
setting  x  =  0  yields 


2  csch  2  C2  = 


r 

-  {1  - 


tanh  C^j' 


From  Equation  (3.34), 


n^  1/2 

tanh  =  (jj— )  coth  C2. 

s 


(3.40) 


Using  this  equation  in  Equation  (3.40)  produces 


2  coth  C2  -  tanh  C2  =  (jjf)  ,  (3.41) 

after  some  rearrangement.  Letting  C2  =  6,  Equation 
(3.41)  can  be  written  as 

2  coth  e  -  tanh  &  *  (j^)  .  (3.42) 

From  Equations  (3.34),  (3.36),  (3.39),  and  (3.42)  the 
normalized  electric  field  is  plotted  in  Figure  3.2  for 
various  substrate-intrinsic  carrier  concentration  ratios. 
On  the  normalized  electric  field  scale  10  is  approximate¬ 
ly  4  x  10^  volts/cm.  The  interface  is  located  at  x  =  0, 

-  X  is  distance  into  the  epitaxial  layer,  and  +  x  is  dis¬ 
tance  into  the  substrate. 

3. 2. 1.3. 2  Surface 

There  is  a  similar  electric  field  at  the  surface. 

The  boundary  condition  for  the  surface  is  that  the  elec¬ 
tron  concentration  is  fixed  at  a  value  n  by  the  surface 
states.  An  expression  for  this  electric  field  similar 
to  Equation  (3.35)  is 

E  (5)  =  -  (i)  ||  .  (3.43) 

Combining  Equations  (3.30)  and  (3.43)  yields  the  electric 


i. 


] 

i 


Figure  3.2  Normalized  electric  field  versus  distance  for 

various  substrate-intrinsic  concentration  ratios, 
N  /n. .  The  interface  is  located  at  x=0,  while 

9  X 

-X  is  the  distance  into  the  epitaxial  layer  and 
+x  is  the  distance  into  the  substrate. 


field  near  the  surface 


E  (5)  =  =1^  [csch  2  +  a)]  (3.44) 

where  a  is  given  by  Equation  (3.29)  and  is  given  by 
Equation  (3.28). 

If  we  combine  Equations  (3.28),  (3.29),  (3.36), 
(3.42),  and  (3.44),  the  electric  field  in  the  epitaxial 
layer  due  to  both  the  surface  states  and  svibstrate  doping 
can  be  obtained.  This  total  electric  field  in  the  epi¬ 
taxial  layer  is  given  by 

(3.45) 


(5)  =  [csch  2  (^  +  a)  +  csch  2  +  6)]  , 


where  a  is  determined  by 


1/2 


tanh  a  =  (j^) 


(3.29) 


6  is  determined  from 

2  coth  6  -tanh  S 


^s 


1/2 


and 


L 


i 


2ee  )cT 
(-^) 
q  n. 


(3.24) 


(3.28) 


With  this  solution  for  the  electric  field,  (5)  and 
N  (C)  can  now  be  determined. 

di 


» 
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3 . 2 . 1 . 4  Impurity  Concentrations 

The  differential  equation  defining  the  impurity  con¬ 
centrations,  Nj  (5) ,  can  be  calculated  from  the  conduction 
current  equations,  the  continuity  equations,  and  the  elec¬ 
tric  field  of  Equation  (3.45).  Combining  Equations  (3.2) 


and  (3.3)  yields 


krfe 


(3.46) 


Changing  Equation  (3.46)  to  the  moving  coordinate  system, 
5=x+vt=x+L,  yields 


g^2  '  Dj  +  kT^  TC  +  kT  j  Dj  3t 


(3.47) 


where  the  upper  -  sign  and  lower  +  sign  are  for  =  +N^ 
and  -N,,  respectively.  Prom  Equation  (3.47)  the  time- 

ci 

independent  differential  equation  for  (C)  is 


j  .  ,  V  3E.  ^  2_ 

^  D”  +  ds  ?  w  a?  ■  °' 


The  solution  of  this  equation  is 


N.  (5) 


(3.48) 


(3.49) 


x/  If  *  il]  «}]. 

n  J 


where  <|»  is  the  electrostatic  potential,  and  (J)^,  E^,  and 
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are  the  potential,  electric  field,  and  donor  or  acceptor 

concentration  at  the  growth  surface.  The  upper  and  lower 

signs  are  for  donors  and  acceptors,  respectively. 

As  an  ex2unple  of  the  results  which  are  obtained  from 

this  model,  we  use  an  epitaxial  growth  temperature  of 

1000  K  {727°C),  where  n^  =  5.8  x  10^®  cm“^  and  »  0.05  ym 

The  ratios  of  the  substrate- intrinsic  carrier  concentra¬ 
tes 

tions,  — ,  and  the  intrrnsic-surface  carrier  concentrations 

ni  "i  2 

are  10  .  The  layer  has  a  thickness  L  =  10  and 

0 

vL  <  <  D j .  These  numbers  correspond  to  the  growth  of  a  thin 
n-layer  on  an  n+  substrate  at  low  growth  rates.  The 
results  are  shown  in  Figure  3.3.  Near  the  growth  surface 
(C  =®  0)  the  electric  field  caused  by  the  surface  states 
tends  to  attract  donor  impurities  and  repel  acceptor  im¬ 
purities  producing  the  impurity  gradients  shown.  Near 
the  layer-substrate  interface  (5  =  L)  the  donor  impurities 
are  repelled  and  the  acceptor  impurities  are  attracted 
by  the  field,  producing  impurity  gradients  in  this  region 
as  well.  The  only  effect  of  changing  layer  thickness 
would  be  to  change  the  size  of  the  units  on  the  horizontal 
axis  (provided  that  vL  remained  less  than  D^). 

Since  donors  in  n-type  GaAs  are  apparently  self- 
compensated  [18],  when  the  sample  is  cooled  to  room  tem¬ 
perature  we  expect  the  net  donor  profiles,  to  be 

pronounced.  Assuming  the  high- temperature  impurity 
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Figure  3.3  Variation  of  donor  and  acceptor  concen¬ 
trations  with  distance  in  the  epitaxial  layer 
at  the  growth  temperature.  The  growth  surface 
is  at  C  =  0  and  the  substrate  is  at  =  L. 
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equllibrium  is  frosen  in,  the  results  of  Figure  3.4  are 
obtained  at  300°  K  for  several  compensation  ratios.  The 
compensation  ratios  indicated  are  values  for  the  middle 
of  the  layer  where  the  gradients  are  zero.  As  can  be 
seen,  the  epitaxial  layer  is  more  heavily  doped  near  the 
growth  surface  and  a  thin  p-type  region  is  obtained  near 
the  layer-substrate  interface. 

In  comparing  these  theoretical  results  with  experi¬ 
mentally  determined  impurity  profiles,  good  quantitative 
agreement  is  not  expected.  This  is  because  in  many  GeJ^s 
epitaxial  growth  situations  the  assmnption  that  vL  is  much 
less  them  is  not  valid  for  donors.  (Little  is  known 
about  the  diffusion  of  the  compensating  acceptors).  Thus, 
a  comparison  with  most  experimental  results  would  require 
the  time  dependent  solution  of  the  problem.  In  this  case 
the  donor  and  acceptor  impurities  are  not  in  equilibrium 
with  the  electric  field  produced  by  the  surface  states  and 
substrate  doping,  and,  qualitatively,  somewhat  different 
behavior  is  expected. 

Near  the  layer-substrate  interface  the  thin  p-region 
should  be  more  pronounced,  since,  when  this  part  of  the 
layer  is  being  grown,  the  fields  from  surface  states  and 
svibstrate  doping  add.  However,  in  either  case  the  model 
produces  a  thin  p-type  region  which  is  in  qualitative 
agreement  with  most  experimental  observations  [7,12,13]. 
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Figure  3.4  Variation  of  the  net  donor  concentration 

with  distance  in  the  epitaxial  layer  at  room 
temperature  for  several  compensation  ratios.  The 
compensation  ratios  indicated  are  values  for 
C  •  0.5L. 
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Also,  in  either  case  the  region  is  expected  to  become 
more  pronounced  as  the  ratio  of  suibstrate  to  layer  con¬ 
centration  increases  [Equation  (3.40)]  in  agreement  with 
experiment  [11] .  In  the  time-dependent  solution,  a  slight 
build-up  of  donors  and  a  depletion  of  acceptors  at  the 
edge  of  the  p- region  in  the  epitaxial  layer  is  expected 
because  of  the  field  direction.  This  would  result  in  the 
type  of  profile  reported  by  Harris  et  al .  [7]. 

Near  the  growth  surface  the  time  dependent  solution 
should  produce  a  heavier  doped  n-region  at  the  surface 
(similar  to  the  time  independent  solution)  with  a  build-up 
of  acceptors  and  a  depletion  of  donors  at  the  edge  of  the 
n-region  in  the  epitaxial  layer.  This  type  of  profile 
near  the  surface  has  not  been  reported  in  the  literature. 
In  differential  capacitance  measurements,  the  surface 
region  where  this  profile  would  be  observed  is  usually 
depleted  by  the  built-in  field  of  the  Schottky  barrier  and 
can  be  measured  only  under  limited  forward  bias .  Although 
impurity  gradients  near  the  surface  have  sometimes  been 
observed  in  forward  bias  measurements,  they  have  generally 
been  considered  as  forward-bias  anomalies  and  have  re¬ 
ceived  little  attention.  If  the  impurity  gradient  model 
is  correct,  surface  impurity  gradients  should  be  common 
in  epitaxial  GaAs  and  could  adversely  affect  device  per¬ 
formance  in  many  applications. 
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The  impurity  gradient  model  discussed  edsove  implies 
that  impurity  gradients  at  the  surface  and  p-type  regions 
at  the  layer-substrate  interface  are  inherent  to  the 
growth  process  when  the  carrier  concentrations  in  the 
substrate,  growing  layer,  and  growth  surface  are  differ¬ 
ent  at  the  growth  temperature.  ThSt  is,  these  problems 
should  be  encountered  in  most  growth  situations  of  prac¬ 
tical  interest.  Considering  only  the  p-type  region  at 
the  interface,  it  is  well  known  that  methods  have  been 
developed  to  eliminate  this  problem.  The  question  is, 
does  the  model  predict  that  these  methods  will  eliminate 
the  problem.  The  most  common  method  of  eliminating  the 
interfacial  region  is  the  use  of  a  buffer  layer  between 
the  substrate  and  the  epitaxial  layer  [8] .  In  practice 
this  constitutes  a  source  of  additional  donors  on  the 
substrate  side  of  the  field,  which  according  to  the  model, 
would  drift  into  the  p-region  and  overcompensate  the 
acceptors.  In  a  similar  manner  the  model  predicts  that 
donors  from  out-diffusion  and  auto-doping  [2]  and 
in  situ  substrate  etching  [7]  would  also  tend  to  compen¬ 
sate  the  p-region,  although  probably  to  a  lesser  extent 
than  the  additional  donors  introduced  during  the  growth 
of  a  buffer  layer.  A  high  arsenic  pressure  during  the 
initial  stages  of  growth  has  also  been  used  to  eliminate 
the  p-region.  The  model  predicts  that  this  could  either 
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eliminate  or  enhance  the  p- region,  depending  upon  the 
specific  nature  of  the  donors  and  acceptors. 

In  conclusion,  the  impurity  gradient  model  dis¬ 
cussed  in  this  section  explains  the  origin  and  prevalence 
of  p-regions  at  n-n+  epitaxial  GaAs  interfaces  and  is 
consistent  with  many  experimental  observations  concerning 
these  regions.  It  also  predicts  an  impurity  gradient  at 
the  outer  surface. 

3.2.2  High  Resistance  Substrate 
3. 2.2.1  Intrinsic  Layer 

From  the  data  in  Section  2.2  chromium-doped  substrates 

are  near  intrinsic  at  1000  K.  If  a  lightly-doped  layer 

which  is  intrinsic  at  the  growth  temperature  is  grown  on 

a  chromium- doped  substrate,  as  can  be  seen  from  Equation 

(3.42),  there  will  be  no  significant  electric  field  at 

^s 

the  interface.  Also^  from  Equation  (3.42),  if  (— )  is 

"i 

approximately  one,  then  3  will  be  infinite,  which  means 
that  the  second  component  of  Equation  (3.45)  will  be  zero. 
Equation  (3.45)  will  then  become 

E(C)  =  -  {  csch  2(j^  +  o)  }  »  (3.50) 

which  is  the  electric  field  due  to  surface  states  only. 

This  electric  field  will  tend  to  attract  donor  impurities 
and  repel  acceptor  impurities  producing  the  Scune  impurity 
gradient  shown  in  Figure  3.4  near  the  surface.  However, 


I 
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there  will  be  no  significant  impurity  gradients  at  the 
interface.  Thus,  for  lightly-doped  epitaxial  layers  only  a 
surface  impurity  gradient  is  expected. 

3. 2. 2. 2  Extrinsic  Layer 

Since  chromium-doped  substrates  are  slightly  p-type 
at  1000  K,  an  extrinsic  layer  at  1000  K  grown  on  a 
chromium-doped  siabstrate  should  have  impurity  gradients 
at  the  interface.  From  Poisson's  equation  the  electron 
concentration  inside  the  epitaxial  layer,  n^,  can  be 
expressed  as 

r  2  2  \ 

+  4n.  } 

n  =  -  , 

®  2 

where  N.  is  the  ionized  donor  concentration,  N  is  the 

CL  cL 

ionized  acceptor  concentration,  and  n^^  is  the  intrinsic 
carrier  concentration.  Of  course,  the  difference  between 
this  electron  concentration  and  the  electron  concentra¬ 
tion  in  the  slightly  p-type  substrate  will  produce  an 
electric  field  at  the  interface.  In  this  electric  field 
donor  impurities  will  be  attracted  towards  the  interface, 
while  acceptors  will  be  repelled.  In  this  manner  an 
impurity  gradient  will  form  to  eliminate  the  electric 
field. 

Because  of  surface  states  the  electron  concentra¬ 
tion  at  the  surface  will  be  lov/er  than  the  electron 
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concentration  inside  the  epitaxial  layer  for  all  but 
heavily-doped  epitaxial  layers.  From  Poisson's  equation 
the  electron  concentration  at  the  surface,  n^,  can  be 
expressed  as 

1/2 

^  *  '"il 

■>0 - ; - 

where  is  the  effective  surface  state  concentration 
sc 

due  to  surface  states.  From  Section  4  this  effective 

surface  state  concentration  is  approximately 

17  -3 

4  X  10  cm  •  With  this  surface  concentration  a  heavily- 

doped  layer  with  a  net  donor  concentration  of 

18  •  3 

4  X  10  cm  will  have  no  significant  electric  field 

at  the  surface.  However,  layers  with  net  donor  concentra- 

17-3 

tions  of  approximately  10  cm  will  have  a  difference 
in  electron  concentration  between  the  surface  and  the 
epitaxial  layer.  This  will  produce  an  electric  field 
which  tends  to  attract  donor  impurities  and  repel  accep¬ 
tor  impurities.  Impurity  gradients  will,  therefore, 
form  such  that  the  concentration  gradients  will  counter¬ 
act  the  electric  field.  This  will  result  in  a  net  donor 

17  -3 

concentration  increase  of  4  x  10  cm  at  the  surface 
to  eliminate  the  electric  field. 


To  summarize,  with  no  outdif fusion  of  impurities 
from  the  substrate,  growth  of  an  extrinsic  layer  on  a 
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chromium- doped  sxabstrate  will  result  in  a  heavier-doped 
region  at  the  interface  similar  to  that  shown  in 
Figure  1.2.  There  will  also  be  a  heavier-doped  region 
located  at  the  surface  due  to  surface  states. 

3 . 2 . 2 . 3  Chromium  Impurity  Gradient 

Tuck  et  al  [72]  have  shown  that  there  is  substantial 
outdiffusion  of  Cr  into  epitaxial  layers  grown  on  Cr- 
loped  substrates.  In  addition,  annealing  experiments  on 
Cr-doped  substrates  with  initially  uniform  chromium  con¬ 
centrations,  have  indicated  a  redistribution  of  atoms 
producing  a  build-up  of  Cr  at  the  surface  [19,21]. 

This  pile-up  of  Cr  at  the  surface  has  also  been  observed 
in  epitaxial  layers  grown  on  C: -doped  substrates  [22,23] 
and  in  ion- implanted  Cr-doped  GaAs  after  annealing  [19,20]. 
This  Cr  gradient  at  the  surface  can  compensate  the  ex¬ 
pected  donor  gradient  calculated  in  Section  3,  thas 
producing  the  zero-bias  capacitance  results  observed  on 
Cr-doped  substrates  (Section  4.2.2). 

In  GaAs,  the  neutral  state  of  chromium  is  Cr^^, 

2+ 

while  the  acceptor  state  is  Cr  [24,73].  Since  both  of 
the  chromium  species  and  the  intentionally-added  donors 
are  positively  charged  they  will  all  drift  in  an  electric 
field  directed  towards  the  surface,  producing  impurity 
gradients.  To  obtain  a  simple  model  for  chromium  in 
GaAs,  we  will  approximate  its  behavior  with  a  single 
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deep  acceptor  near  midgap  [74,75].  It  is  well-known, 
however,  that  its  effects  are  somewhat  more  complicated 
than  this  model  [73,76]. 

To  excuaine  the  formation  of  a  chromium  impurity  gra¬ 
dient  at  the  surface,  the  following  assumptions  ire  made: 

(1)  The  carrier  concentration  of  the  growing  layer, 
n^^,  is  determined  by  the  growth  temperature  and  not  by 
the  ionized  donors,  N, ,  and  acceptors,  N  ,  incorporated 

Cl  A 

during  the  growth  process  or  the  chromium  concentration, 
^Cr'  O'itdif fusing  from  the  substrate.  That  is,  n^  is 
greater  than  N^,  and 

(2)  The  carrier  concentration  at  the  growth  surface, 
n^,is  fixed  by  the  surface  states  and  is  less  than  n^. 

(3)  The  product  of  growth  rate,  v,  and  epitaxial 
layer  thickness,  L,  is  much  less  than  the  impurity  dif¬ 
fusion  coefficients,  D ^ ,  where  j  represents  donors,  ac¬ 
ceptors,  or  chromium. 

These  assumptions  are  similar  to  those  used  before, 
so  that  the  resulting  electric  field  near  the  surface  is 

E(5)  =  -  csch  2(^  +  a),  (3.51) 

where  a  is  determined  by 

tanh  a  =  . 


(3.52) 
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The  spatial  variable  C  is  obtained  from  the  moving  coor¬ 
dinate  system 

5=x+vt=x+L,  (3.53) 

where  C  is  the  distance  from  the  growth  surface,  and 
-X  is  the  distance  from  the  substrate.  From  this  elec¬ 
tric  field  the  electrostatic  potential,  $,  is 


<()  (5)  =  ^  In  1  tanh  (i  +  a)  |  .  (3.54) 

To  calculate  the  behavior  of  the  donor,  acceptor, 
and  chromium  atoms  produced  by  this  electric  field,  we 
use  the  current  equations  for  ionized  donors,  acceptors, 
and  both  ionized  and  neutral  chromium  atoms. 


Jj  =  Zqy^NjE  -  ZqDj  ; 


(3.55) 


and  the  continuity  equations; 


3N.  ,  3J. 

Tt  “  +  ^  ^ 

The  upper  -  signs  in  Equations  (3.55)  and  (3.56)  are  for 
the  positively  charged  donors  and  chromium  atoms,  while 
the  lower  +  signs  are  for  the  negatively  charged  accep¬ 
tors.  The  value  of  Z  in  these  equations  is  one  for  the 
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singly  ionized  donors  and  acceptors,  and  either  two  or 
three  for  the  chromium  atoms. 

Using  Equation  (3.55)  in  Equation  (3.56)  and 
Einstein's  relation. 


‘j  “  It  °j  ' 


(3.57) 


yields 


1  M  .  ’S  ,  V  ^  j  za  » 

Dj  ft  “  D"  ^  W  TC  * 


(3.58) 


Under  the  assumption,  vL  <<  D ^ ,  the  system  is  in  diffusive 
equilibrium  and  the  time- independent  differential  equa¬ 
tion  for  Nj  (C)  is 


d^N^ 


d? 


(3.59) 


The  solution  of  Equation  (3.59)  is 


.  exp  +  Z2_<ii) 
N  'D.  kT  ' 

0  D 


f  .,*'’♦0,  ,r  V  T  .  1  “"d 

|exp  -<  XT  '  L  Dj  XT  *  05  J 

X  /  exp  (-  gi  t  ^a|i5i)  dc}  , 


(3.60) 


where  the  electrostatic  potential  4)  (5)  is  given  by 


I 

» 
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Equation  (3.54)  and  are  the  potential, 

electric  field,  and  donor,  acceptor  or  chromium  con¬ 
centration  at  the  growth  surface.  The  upper  signs  are 
for  donors  and  chromium  atoms  while  the  lower  signs  are 
for  acceptors.  The  value  of  Z  is  1  for  ionized  donors 
and  acceptors,  2  for  ionized  chromium  atoms,  and  3  for 
neutral  chromium  atoms. 

Quantitatively,  we  consider  an  epitaxial  growth 

temperature  of  1000  K  (727°C) .  At  this  temperature 
16 

nj^  =  5.8  X  10  cm  ,  and  =  0.05  pm.  For  purposes 

of  illustration  we  assume  n./n^  =  10  and  L/,  =  10. 

i  0 

The  results  of  Equation  (3.60)  for  these  conditions  are 
shown  in  Figure  3.5.  Near  the  growth  surface  (5*0)  the 
electric  field  caused  by  the  surface  states  tends  to 
attract  donor  and  chromium  impurities  and  repel  acceptor 
impurities,  producing  the  impurity  gradients  shown. 

Assuming  that  the  high-temperature  impurity  profiles 
are  frozen  in  when  the  growth  process  is  terminated 
and  the  epitaxial  layer  is  cooled  to  room  temperature, 
the  net  donor  profile,  ^  shown  in 

Figure  3.6  for  various  concentrations  of  chromium  at 
5  =  5Lj^  in  the  epitaxial  layer.  The  compensation  ratio 
(N^  +  N^^) /Nj  was  fixed  at  0.25  for  various  N^^  by 
adjusting  N  .  This  net  donor  profile  shows  the  forma- 
tion  of  a  high  resistivity  or  p-type  region  at  the 


Figure  3.6  Variation  of  th( 
at  room  tempera' 


trations .  The 
ratios  are  val 
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svirface  if  there  is  chromium  in  the  epitaxial  layer. 

If  the  chromium  concentration  is  only  3  per  cent  of  the 
donor  concentration  at  C  ^  a  region  of  negative  net 

donor  concentration  will  result  at  the  svirface  (assuming 
that  all  of  the  chromium  atoms  are  electrically  active) . 

It  has  been  found  that  not  all  of  the  chromium  in 
GaAs  is  electrically  active  as  a  deep  acceptor.  Part  of 
the  total  chromium  atoms  may  be  either  neutral  or  at 
interstitial  sites,  acting  as  donors  [77].  In  this  man¬ 
ner  it  may  compensate  itself  [78].  Another  possibility 
is  that  the  high  resistivity  seen  in  Cr-doped  GaiAs  is 
indirectly  caused  by  the  chromium  [76].  That  is,  the 
deep  acceptor  may  be  a  Cr-0  complex  [76]. 

In  conclusion,  the  lower  zero-bias  capacitance  ob¬ 
served  for  the  as-grown  portion  versus  the  etched  portion 
of  an  epitaxial  layer  grown  on  Cr-doped  substrates,  can 
be  explained  by  a  chromium  impurity  gradient  formed  due 
to  surface  states.  This  chromium  impurity  gradient  can 
compensate  the  donor  impurity  gradient  and  the  donor 
concentration  of  the  epitaxial  layer,  causing  the  de¬ 
pletion  width  under  the  Schottky  barrier  to  increase  to 
uncover  the  same  amount  of  charge.  Since  capacitance 
is  inversely  proportional  to  the  depletion  width,  the 
capacitance  will  decrease. 

This  formation  of  a  chromium  impurity  gradient  due 


to  surface  states  can  also  explain  the  observed 
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redistribution  of  chromium  f 19, 20, 21]  when  Cr-doped 

GeiAs  is  annealed.  The  observed  pile-up  at  the  surface 

0 

seen  by  Evans  et  al.  [20]  is  typically  under  1000  A, 
which  is  the  order-of-magnitude  predicted  by  the  model. 
It  can  also  explain  the  abnormally  high  diffusion  of 
chromium  towards  the  surface  during  the  annealing  of 
Cr- implanted  G2As  [79].  Thus,  the  chromium  gradients 
observed  in  these  annealing  experiments  could  be  caused 
by  siirface  states  rather  than  strain  from  the  encapsu- 
lant  as  has  been  suggested  [21 ] . 


-83- 


4.  SURFACE  STATE  EFFECTS 

4.1  SURFACE  IMPURITY  GRADIENTS 

4.1.1  Analysis 

The  model  presented  in  Section  3.2  predicts  that 
there  will  be  impurity  gradients  at  the  surface  of  an 
epitaxial  layer  due  to  surface  states.  These  impurity 
gradients  result  in  the  net  donor  profile  shown  in 
Figure  4.1  for  an  as-grown  epitaxial  layer.  Unfortunate¬ 
ly,  it  is  difficult  to  obtain  a  direct  measurement  of 
this  profile  near  the  surface.  With  differential  capaci¬ 
tance  measurements  this  profile  is  already  depleted  at 
zero-bias  by  the  built-in  field  of  the  Schottky  barrier. 
Although  in  principle  it  is  possible  to  probe  this  region 
with  foward-bias  measurements,  the  much  larger  current 
and  resulting  lower  Q  of  the  circuit  produce  questionable 
results.  This  prevents  a  direct  measurement  of  the 

1  g  —3 

net  donor  gradient  for  all  doping  levels  below  7  x  10  cm 
However,  an  indirect  measurement  of  this  net  donor  pro¬ 
file  can  be  obtained  from  the  zero-bias  capacitance  of 
the  Schottky  barrier.  When  a  Schottky  barrier  is  formed 
on  the  surface  of  an  as-grown  epitaxial  layer  with  a  sur¬ 
face  gradient,  a  certain  quantity  of  charge  under  the 
barrier  will  be  uncovered  due  to  surface  states.  The 
measured  zero-bias  capacitance  corresponds  to  a  deple¬ 
tion  width  Wq^.  If  the  surface  region  of  the  GeiAs  with 


this  net  donor  gradient  is  then  etched  away,  so  that  the 
net  donor  profile  is  uniform,  the  depletion  width  under 
the  Schottky  barrier  must  increase  to  to  uncover  the 
scuae  total  charge.  The  excess  charge  produced  by  the 
impurity  gradient  is  the  area  under  the  as-grown  profile 
in  Figure  4.1  out  to  the  as-grown  depletion  width, 
minus  the  area  under  the  etched  profile  also  out  to  Wq^ 
Assiiming  the  surface  state  densities  are  the  Scime,  this 
equals  the  area  increase  under  the  etched  profile  in 
Figure  4.1  caused  by  the  increase  of  the  zero-bias 
depletion  width  for  the  etched  epitaxial  layer  from 
to  Wqj..  This  increase  in  depletion  width,  of  course, 
corresponds  to  a  decrease  in  zero-bias  capacitance  for  the 
etched  layer. 

For  an  as-grown  epitaxial  layer  with  impurity  gra¬ 
dients  at  the  surface,  the  charge  uncovered  beneath  a 
Schottky  barrier  is 

{n^  (5)  -  (5)}  dS,  (4.1) 

0 

where  (5)  and  (5)  are  taken  from  Section  3.2.1  as 

Nd  (5)  =  (5L^)  coth^  (^*2) 

(C)  =  (5L. )  tanh^  +  a) . 

a  ax  Xj» 
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In  these  equations 


2££kT  and 

q  n. 


-1 

n  =  tanh  ^  (^) 
“i 


Using  Equation  (4.2)  in  Equation  (4.1)  yields  a  total 
uncovered  charge  given  by 


Qqa  = 


j. 

I 


,  (5L.  )  coth^  (^  +  a) 
a  1 


(4.3) 


-N^(5L^)  tanh^ 


where  is  calculated  from 


W, 


ee^A 


OA  C 


OA 


(4.4) 


Cq^  is  defined  as  the  as-grown  zero-bias  capacitance. 
Performing  the  integration  in  Equation  (4.3)  yields 


OOA  -  '^“■1  {  < W  ^  <<•  =  > 

w 

+  ^  coth  a  -  coth  (-^^  +a)j 

w 

-  ^tanh  a  -  tanh  ^  ' 

where  =  N,  (5L. )  and  N,  =  N,  (5L. ) . 

Q  Q  X  a  »  X 
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For  the  epitaxial  layer  with  the  impurity  gradients 
etched  away,  the  uncovered  charge  is 


QqE  =  WqE'  (4.6) 

where  is  obtained  from  the  etched  zero-bias  capaci¬ 
tance  ,  Cq^ ,  as  in  Equation  (4.4).  Since  the  uncovered 
charges,  and  Qq£/  are  equal,  then  the  right-hand  sides 
of  Equations  (4.5)  and  (4.6)  are  also  equal.  Combining 
Equations  (4.5)  and  (4.6)  gives 


qALi  ^  =  qALj^  { 

+  ^  coth  a  -  coth  (-^^  +  a)  j 

-  ^  tanh  a  -  tanh  (-^^  ^  ]  1  * 


^OA 

Experimentally,  -j—  is  relatively  large  and 
coth  +  a)  =  tanh  (-^^  +  a)  =  1. 


Equation  4.7  becomes 


(4.7) 


W, 


OE 


W 


OA 

^i 


+  1)  = 


N,  coth  a  -  N  tanh  a 

d  dL. 

d  A 


Also,  N  tanh  a  is  much  smaller  than  N,  coth  a. 
a  d 


(4.8) 


5' 
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Thus,  Equation  (4.8)  is  given  approximately  as 


N 

coth  a  =  (— 
Solving  Equation 


d-\.  f  ”oE 
“d  ^  ^i  " 

(4.9)  for  a  yields 


(4.9) 


(4.10) 


Thus,  given  the  zero-bias  capacitances  and  assuming  a 
reasonable  compensation  ratio  the  surface  concentra¬ 

tion  n^  needed  to  produce  the  impurity  gradients  can  be 
calculated  from  Equation  (4.10)  and  the  boundary  condition 


n^  =  n^  tanh^  a.  (4.11) 

With  this  surface  electron  concentration  and  Poisson's 
equation , 


VE  =  Z  .  N.  (4.12) 

where  N.  =  +  N,,  -  N,,  +  p,  -  n, 

J  Q  a 

an  effective  surface  state  concentration  N  which  equals 

SC 

ee 

VEq  can  be  calculated.  Finally,  the  surface  state 
density,  N  ,  needed  to  create  this  surface  concentration, 

S 


can  be  obtained. 
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4.1.2  Experiment 
4 . 1 . 2 . 1  Heavily-Doped  Substrate 

Epitaxial  layers  were  grown  on  n'*’  GeiAs  siibstrates 
in  an  AsCl2~Ga-H2  vapor-phase  reactor.  The  substrates 
had  a  crystallographic  orientation  2°  off  <100>  to  pre¬ 
vent  the  formation  of  hillocks  [80].  The  doping  was 
controlled  by  the  addition  of  H2S  to  the  reactor.  In 
this  manner,  the  epitaxial  layers  were  grown  to  a  thick¬ 
ness  of  between  4  and  7  ym  with  doping  concentrations  in 
the  10^^-10^®cm~^  range.  The  samples  were  then  removed 
from  the  reactor  and  a  portion  of  each  wafer  was  taOcen 
and  etched  in  either  H2SO^ :H202 :H20  (5:1:1)  or  hot  Auro- 
strip  such  that  2000-3000  R.  of  GaAs  was  removed.  Sample  4 
of  Table  4 . 1  was  the  only  sample  to  be  etched  in  Aurostrip 
while  Samples  1,  2,  and  3  were  etched  in  H2S0^ :H202 SH2O. 
Gold  Schottky  barriers  were  formed  by  plating  through 
a  photoresist  mask  (Scunple  3)  ,  evaporating  (Scunple  2)  , 
and  sputtering  (Scunples  1  and  4) .  The  values  of  the  zero- 
bias  capacitance  measured  for  both  the  as-grown  and  etched 
portions  of  Samples  1  through  4  are  given  in  Table  4.1. 
These  capacitances  are  normalized  to  unit  area.  Typically, 
15-25  diodes  were  measured  on  both  sections  of  the  wafers 
with  the  average  values  and  standard  deviations  shown  in 
Table  4.1.  Reverse  bias  differential  capacitance  measure¬ 
ments  on  each  section  indicated  uniform  carrier 
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Table  4.1  Zero-bias  capacitance  for  etched  and 

as-grown  on  n^  substrates  and  the  surface 
state  densities  at  1000  K  needed  to 
produce  these  data. 


Scunple 

Zero-bias  Capacitance 

Surface  State  Density 

Cq  (nf/cm^) 

(cm"^) 

s 

As-Grown 

Etched 

At  1000  K 

1 

16.7±0.4 

13.7±0.3 

5.3x10^^ 

2 

17.3±0.5 

14.6±0.4 

4.8x10^^ 

3 

19.1±0.9 

15.4±0.4 

5.6x10^^ 

4 

91.3±2.0 

74.012.5 

4.5x10^^ 

concentrations  beyond  the  surface  regions.  Using  Equa¬ 
tions  (4.10)  and  (4.11)  and  assuming  a  compensation  ratio, 
N^/N^,  of  0.30,  the  surface  concentrations  were  calculated 
for  the  seunples  listed  in  Table  4.1.  Then,  using  Equa¬ 
tion  (4.12),  the  effective  surface  state  concentrations 
were  obtained.  The  surface  state  densities  at  1000  K  for 
Seunples  1  though  4  are  also  listed  in  Table  4.1.  From 
these  data  the  average  surface  state  density  for  GaAs 
at  1000  K  was  5.1  x  lO^^cm"^. 

4 . 1 . 2 . 2  Chromium- Doped  Substrate 

Epitaxial  layers  were  grown  on  Cr-doped  GaAs  substrates 
in  an  AsCl2-Ga-H2  vapor-phase  reactor  [39,40].  The  sub¬ 
strates  had  a  crystallographic  orientation  2°  off  <100>. 

14  -3 

With  a  background  concentration  around  10  cm  ,  the 
doping  of  the  epitaxial  layers  was  controlled  by  the  addi¬ 
tion  of  tin  to  the  gallium  melt  [81] .  In  this  manner, 
the  epitaxial  layers  were  grown  to  a  thickness  of  10  and 
20  wm  with  doping  concentrations  in  the  lO^^-lO^®  range. 

As  before,  the  scunples  were  removed  from  the  reactor  and 

a  portion  of  each  was  etched  in  H2SO^ ;H202 :H20  (5:1:1) 

0 

such  that  2000-3000  A  of  GaAs  was  removed.  Of  the  seun¬ 
ples  listed  in  Table  4.2,  gold  Schottky  barriers  were 
formed  by  plating  for  Samples  1  and  2  and  by  evaporation 
for  Seunples  3,  4,  and  5.  The  zero-bias  capicitance  values 
normalized  to  unit  area,  are  listed  in  Table  4.2.  These 
data  reflect  impurity  gradients  different  from 
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those  obtained  for  layers  grown  on  heavily-doped  substrates. 
That  is,  the  as-grown  zero-bias  capacitance  values  are 
lower  than  the  etched  capacitance  values,  indicating  a 
higher  resistivity  region  at  the  surface. 

One  possible  explanation  could  be  the  creation  of 
arsenic  vacancies  near  the  surface  during  furnace  cool 
down.  However,  arsenic  vacancies  are  believed  to  be  do¬ 
nors  [82,83,84]  which  would  cause  an  n+  region  at  the  sur¬ 
face  rather  tham  the  observed  higher  resistivity  region. 

Because  of  the  high  diffusion  rate  of  chromiuin  from 
the  substrate,  a  more  likely  explanation  of  the  data  in 
Table  4.2  is  a  pile-up  of  chromium  at  the  surface  due  to 
the  electric  field  produced  by  surface  states.  Since 
chromium  in  GaAs  is  positively  charged  [24,73],  the 
chromium  should  drift  towards  the  surface,  compensating 
the  donors  also  drifting  towards  the  surface. 

As  was  demonstrated  in  Section  3. 2. 2. 3,  this  double 
drift  process  under  certain  conditions  can  produce  a 
high- resistivity  region  near  the  surface.  The  compen¬ 
sation  of  the  donor  impurity  gradient  and  the  donor  con¬ 
centration  of  the  epitaxial  layer  by  the  chromium 
gradient  will  cause  the  depletion  width  under  a  Schottky 
barrier  to  increase  to  uncover  the  seune  eunount  of  charge. 
Since  capacitance  is  inversely  proportional  to  the  de¬ 
pletion  width,  the  zero-bias  capacitance  will  decrease, 
producing  results  similar  to  the  data  in  Table  4.2. 
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4.2  SELF-COMPENSATION  EFFECT 

When  the  donor  and  acceptor  concentrations  are  far 
below  the  Intrinsic  carrier  concentration,  the  impurities 
are  screened  from  each  other  by  the  intrinsic  carrier 
concentration  during  incorporation  into  the  epitaxial 
layer.  This  produces  a  constant  compensation  ratio, 

N^/N^,  with  increasing  total  impurity  concentration, 

N^+N^.  However,  when  the  impurity  concentrations,  N^  and 
N_,  reach  the  Scune  order-of-magnitude  as  the  intrinsic 
carrier  concentration,  the  impurities  are  no  longer  ef¬ 
fectively  screened.  The  interaction  between  the  impuri¬ 
ties  then  affects  the  incorporation  of  the  donors  and 
acceptors  causing  the  compensation  ratio  N^/N^  to  in¬ 
crease.  At  high  total  impurity  concentrations  the 
acceptor  concentration  will  equal  the  donor  concentration. 
4.2.1  Experiment 

Figure  4 . 2  shows  the  donor  and  acceptor  concentra¬ 
tions  versus  the  total  impurity  concentration  for  a 
number  of  epitaxial  layers  doped  with  column  IV  impuri¬ 
ties  [18].  These  data  exhibit  the  expected  increasing 
compensation  ratio  at  high  total  impurity  concentrations. 

These  epitaxial  layers  were  grown  on  Cr-doped  semi- 
insulating  GaAs  substrates  in  an  AsCl^-Ga-H2  vapor  phase 
reactor  [39,40].  The  substrates  had  a  crystallographic 
orientation  of  <211>  Ga  because  this  orientation  in¬ 
corporates  the  smallest  number  of  residual  impurities  [42] . 


OR  (cm"®) 


Nj+Ng  (cm"®) 


Figure  4.2  Variation  of  the  donor  concentration  (N^)  and  the 
acceptor  concentration  (N  )  with  the  total 

cl 

ionized  impurity  concentration  (N,+N  )  for  the 

Cl  ^ 

column-IV  impurities:  Si,  Ge,  and  Sn.  The  lines 

11  -2 

are  calculated  taking  into  account  6.3  x  10  cm 
surface  states. 


the 
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With  a  backgrotind  concentration  around  10  cm  , 

doping  was  controlled  by  the  addition  of  Si,  Sn  or  Ge 

to  the  gallium  melt  [85].  In  this  m2mner  epitaxial  layers 

14 

were  grown  with  carrier  concentrations  between  5  x  10 
and  5  x  lO^^cm  The  donor  concentrations,  N^,  and  the 
compensating  acceptor  concentrations,  N_,  were  determined 

cl 

from  analyses  of  Hall  and  resistivity  measurements  [45] . 

As  can  be  seen  in  Figure  4.2,  both  and  in¬ 
crease  linearly  with  respect  to  N .  +  N  ,  (N^/N,  remains 
constant)  at  low  total  impurity  concentrations  as  ex¬ 
pected.  At  high  concentrations,  however,  the  compensa¬ 
tion  ratio  does  not  begin  increasing  until 

is  about  an  order-of-magnitude  larger  than  the  intrinsic 
concentration.  Therefore,  something  with  a  higher  con¬ 
centration  than  the  intrinsic  carriers  is  effectively 
screening  the  donor-acceptor  interaction. 

Since  surface  states  increase  the  hole  concentration 
in  n-type  material  at  the  growth  surface  [86] ,  surface 
states  could  effectively  screen  the  impurities  to  con¬ 
centrations  higher  than  the  intrinsic  carrier  concentra¬ 
tion.  Casey  et  al.  [87]  and  Zechauer  and  Vogel  [88]  have 
used  surface  band  bending  [86]  to  account  for  several 
previously  unexplained  impurity  incorporation  effects 
in  the  liquid-phase  epitaxial  growth  of  GaAs,  From 
Section  4.1.2,  we  see  that  band  bending  due  to  surface 
states  is  applicable  to  vapor  phase  epitaxy  as  well. 
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I£  the  screening  effect  of  surface  states  Is  taken  into 
account,  the  solid  lines  in  Figure  4.2  can  be  fit  to 
the  data.  This  will  be  discussed  in  detail  in  Section 
4. 2. 2. 3. 

4.2.2  Analysis 

4. 2. 2.1  Mass  Action  Laws 

To  model  the  incorporation  of  impurities  into 
epitaxial  GaAs,  mass  action  laws  for  dilute  solutions  [89] 
are  used.  The  first  quasi-chemical  reaction, 

0  =  e  +  h"*" , 

is  for  the  generation  of  holes  and  electrons.  The 
associated  mass  action  law  is 


where  n^  is  the  equilibrium  constant  which  is  strongly 
temperature-dependent.  Assuming  Schottky  disorder,  the 
quasi-chemical  reaction  and  associated  mass  action  law 
for  the  formation  of  stoichiometric  defects  are 
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where  P^g  is  the  partial  pressure  of  as^  in  the  reactor. 

The  addition  of  an  impurity,  L,  into  the  epitaxial 
layer  is  represented  by 


and 


L(s)  +  Vj  =  Lj 


[L^l  =  [Vjl, 


(4.16) 


where  the  j  may  be  either  a  gallium  or  arsenic  site  and 
Nt  is  the  total  concentration  of  the  impurity.  Assuminc 

id 

a  donor  impurity,  the  ionization  reaction  is  given  by 


and 


L.  =  L.  .  e 


nN^  =  [L.]. 


(4.17) 


Considering  the  self-compensation  of  doped  n-type 
GaAs  [18,90],  there  is  now  substantial  evidence  that  the 
compensating  acceptors  [18,83,91-98]  are  impurity-defect 
complexes.  Because  of  the  linear  increase  of  with 
total  added  donor  atoms  for  the  data  in  Figure  4.2, 
this  complex  must  have  a  donor  atom  associated  with  it. 
We  assume  that  the  complex  is  a  donor  atom  associated 
with  a  vacancy  on  the  same  sub-lattice  as  the  donor  atom. 


-99- 


This  acceptor  complex  is  represented  by 
L(s)  +  2  Vj  =  (Lj  Vj) 

and  (4.18) 

[(Lj  V.)]  =  K2N^  [V.]2. 

The  ionization  of  this  acceptor  complex  is  given  by 

(LjVj)  =  (LjVj)"  +  h"^, 

and  (4.19) 

p  “a  -  K,  [(L.Vj)l. 

It  should  be  pointed  out,  however,  that  for  constant 
temperature  and  partial  pressure  of  arsenic,  the  analysis 
in  Sections  4. 2. 2. 2  and  4. 2. 2. 3  is  independent  of  the 
type  of  self-compensating  acceptor  assumed.  An  amphoteric 
acceptor  [17] ,  or  another  donor-vacancy  complex  [94-98] 
will  produce  the  same  results. 

Using  Equations  (4.14,  4.15,  4.16,  and  4.17),  the 
ionized  donor  concentration,  N^,  is 


Nd  ={ 


^d^l^As 


n  p 


1/4 


As 


(4.20) 


for  incorporation  of  the  donors  on  the  arsenic  sub¬ 
lattice.  For  incorporation  of  the  donors  on  the  gallium 
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sub- lattice  the  ionized  donor  concentration  is  given  by 


d  L  n 


Pas 


]  “l  • 


(4.21) 


For  constant  temperature  and  arsenic  partial  pressure 
both  Equations  (4.20)  and  (4.21)  become 


^d  =  IT  ^  N^) 


(4.22) 


where  the  total  impurity  concentration, 

In  a  similar  manner,  the  acceptor  concentration, 

N_,  can  be  calculated  using  Equations  (4.13,  4.14,  4.15, 

cL 

4.18,  and  4.19)  resulting  in  the  following: 


n  K,K-K,„2 
N  =  a  2  As 

“a  2  1/2  ‘L' 

■^i  Pas 


(4.23) 


for  complexes  on  the  arsenic  sub-lattice,  and 


”a  n  2  ^2  L' 

1  As 


(4.24) 


for  complexes  on  the  gallium  sub-lattice.  For  constant 
arsenic  partial  pressure  and  temperature.  Equations  (4.23) 
and  (4.24)  can  be  represented  by 


Na  =  "  (^d 

where  N ,  +  N  *  . 

a  a  Xi 


(4.25) 
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The  net  donor  concentration  can  then  be  expressed 
as 


K  **  K  n 

'Va>  -{-i - 

"i 

using  Equations  (4.25)  and  (4.22).  When  the  impurity  con¬ 
centrations  are  small  compared  to  n^  the  layer  is  intrinsic 
and  Equation  (4.26)  becomes  »  which  pro¬ 

duces  constant  compensation.  As  the  carrier  concentration 
increases  at  high  doping  levels,  K^/n  decreases  while 
K^n/n.  increases  causing  the  compensation  ratio  to  increase. 

cL  X 

4. 2. 2. 2  Charge  Neutrality 

If  we  assume  no  surface  states,  then  there  is  no 
electric  field  at  the  surface  and  charge  neutrality 
must  be  maintained.  Thus,  the  neutrality  condition  is 


n  +  N  =  p  +  N 

a  ^  r 


(4.27) 


for  the  epitaxial  layer.  This  equation  can  be  rewritten  as 
_  *  4n2,  1/2 

n  =  - ,  (4.28) 


by  using  Equation  (4.13).  From  this  equation,  if 

N  <<  n.  then  the  electron  concentration,  n,  is  equal 

Cl  a  X 

to  the  intrinsic  carrier  concentration. 

16  —  3 

For  a  total  impurity  concentration  below  10  cm  , 
the  carrier  concentration  is  within  5  percent  of  the 
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intrinsic  concentration.  From  the  data  in  Figure  4.2, 
the  equilibrium  constant  in  Equation  (4.22)  can  be 
calculated.  The  calculated  value  of  Kj^  is  4.90  x  lO^^cm”^ 
with  a  standard  deviation  of  0.14  x  lO^^cm”^  for  29  data 
points. 

Similarly,  from  the  data  of  Figure  4.2,  the  equilibri¬ 
um  constant  of  Equation  (4.25)  is  calculated  to  be 

Ci 

16  *“1 

1.10  X  10  cm  .  The  standard  deviation  of  this  value 
is  0.82  X  lO^^cm”^  for  29  data  points.  From  the  data  in 
Figure  4.2,  the  calculated  equilibrium  constants,  and 
K',  and  Equation  (4.28),  the  net  donor  concentration  can 
be  calculated  from  Equation  (4.26)  for  no  surface  states. 
This  calculated  net  donor  concentration  is  the  dashed 
line  shown  in  Figure  4.3.  As  can  be  seen,  the  experimen¬ 
tal  data  do  not  agree  with  this  calculated  curve.  There¬ 
fore,  some  charged  species  at  a  higher  concentration  than 
the  intrinsic  carrier  must  be  screening  the  donor  and 
acceptor  impurities  as  they  are  being  incorporated  into 
the  epitaxial  layer. 

Hurle  [99]  has  proposed  a  model  to  explain  this  be¬ 
havior.  His  model  is  based  on  a  large  concentration  of 
arsenic  vacancies  due  to  Frenkel  disorder  on  the  arsenic 
sub-lattice.  For  his  model  the  arsenic  vacancy  is  taken 
to  be  a  donor.  Thus,  at  high  temperatures  the  electron 
concentration,  n,  equals  the  ionized  arsenic  vacancy 


(cm**) 


Figure  4.3  Net  donor  concentration  versus  total 

ionized  impurity  concentration  (N^+N^) .  The 

dashed  line  is  calculated  for  zero  surface 

states.  The  solid  line  is  calculated  for 
11  —2 

6.3  X  10  cm  surface  states.  The  layers 
are  doped  with  the  column- IV  impurities: 

Si  (•)  ,  Ge  (▲)  ,  and  Sn  (■)  . 
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concentration  [V^g]  which  is  much  greater  than  the  intrinsic 
carrier  concentration,  n^.  This  concentration  of  electrons 
and  arsenic  vacancies  screen  the  donors  and  acceptors  as 
they  are  incorporated  into  the  epitaxial  layer.  At  1000  K 
this  model  requires  an  electron  concentration  of  6  x  10^^ 
cm”^  due  to  the  arsenic  vacancies  to  explain  the  experi¬ 
mental  data.  However,  the  measured  electron  concentration 
at  1000  K  from  the  high  temperature  Hall  measurements  in 
Section  2  is  an  intrinsic  concentration  of  5.8  x  lO^^cm”^. 
Thus,  the  basic  assumption  of  Hurle's  model  is  incorrect. 

4. 2. 2. 3  Surface  States 

The  effect  of  surface  states  is  another  explanation 
for  the  screening  of  impurities  above  the  intrinsic  carrier 
concentration.  As  shown  in  Section  3  these  surface  states 
will  produce  an  electric  field  at  the  surface  of  the  epi¬ 
taxial  layer.  Because  of  this  electric  field  Poisson's 
equation, 


-^VE  =  +  p  -  -  n. 


(4.29) 


must  be  used  instead  of  charge  neutrality.  The  gradient 
of  the  electric  field  at  the  surface  corresponds  to  a  sur¬ 
face  concentration,  N  ,  due  to  the  surface  states.  Thus, 

S  w 

at  the  surface  Equation  4.29  can  also  be  written  as. 


N  +N  +n=N,  +p. 
sc  a  d 


(4.30) 


t 
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If  the  interaction  of  the  impurities  with  the  field  is 

assumed  to  be  at  the  surface  only,  then  N  can  be  taken 

sc 

as  a  constant.  Equation  4.30  can  solved  for  the  elec¬ 
tron  concentration  as. 


‘"d  -  "a  -  "sc’  “•'d  -  “’a  - 


(4.31) 


For  impurity  concentrations,  and  N^,  less  than  the  in¬ 
trinsic  concentration.  Equation  4.31  becomes. 


n 


4n2]V2 


2 


(4.32) 


which  is  constant  for  constant  temperature. 

To  solve  for  the  equilibrium  constants  and  K^, 

the  data  in  Figure  4.2  below  ■*'  “  lO^^cm  ^  and  an 

assumed  N  are  used.  With  Equations  4.22,  4.25,  4.31 
sc 

and  these  equilibrium  constants,  the  N  needed  to  fit  the 

5  C 

data  above  N,  +  N  =  lO^^cm”^  is  determined.  By  iteration 
d  a 

all  the  data  in  Figure  4.2  are  then  fit  with  a  self-con¬ 
sistent  surface  state  concentration. 

In  this  manner,  a  surface  state  concentration  of 
5.0  X  lO^^cm”^  was  obtained.  The  calculated  equilibrium 


constants  were 
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»  5.58  X  10^^  +  .08  X  lO^^cm"^ 

and 

K'  =  9.62  X  10^®  +  .78  X  lO^^cm"^, 

where  the  standard  deviations  are  for  29  data  points. 

With  these  constants.  Equations  (4.26)  and  (4.31)  and  a 

surface  state  concentration  of  5.0  x  10 the  solid 

line  plotted  in  Figure  4.3  was  determined.  The  surface 

state  density  to  produce  this  surface  concentration  was 
11  -2 

6.3  X  lO-^-^cm 

4.3  CONCLUSIONS 

From  the  previous  two  sections  it  is  apparent  that 
surface  states  can  effect  the  incorporation  and  distribu¬ 
tion  of  impurities  in  vapor-phase  epitaxy.  Prom  Section 

4. 1.2.1  the  average  surface  state  density  at  lOOOK  needed 

to  produce  the  surface  impurity  gradients  observed  was 
11  -2 

5.1  X  10  cm  .  This  surface  state  density  corresponds  to 
a  barrier  height  of  0.17eV  for  the  vapor-solid  interface 
at  lOOOK. 

From  Section  4. 2. 2. 3  the  surface  state  density  neces- 

11  -2 

sary  to  screen  the  impurities  was  6.3x10  cm  .  A  bar¬ 
rier  height  of  0.20eV  is  produced  by  this  surface  state 
density.  This  surprisingly  good  agreement  between  two 
completely  different  experiments  tends  to  confirm  the 


effects  of  surface  states  on  the  growth  of  epitaxial  layers. 

The  electric  field  produced  by  these  surface  states 
may  also  explain  why  G2J^s  epitaxial  layers  are  self-compen¬ 
sated.  That  is  the  field  could  be  expected  to  retard  the 
formation  of  donors  from  donor  atoms  and  enhance  the  forma¬ 
tion  of  acceptors  from  donor-atom-vacemcy  complexes,  thus 
compensating  the  layers. 


-108- 


5.  SUMMARY  AND  CONCLUSIONS 

The  impurity  gradient  model  presented  in  this  report 
shows  that  impurity  gradients  are  inherent  to  the  growth 
process  of  epitaxial  GaAs.  That  is,  if  all  other  pos¬ 
sible  mechanisms  for  the  fojnnation  of  impurity  gradients 
were  eliminated,  in  most  epitaxial  growth  situations  of 
practical  interest  these  regions  would  still  be  observed. 
This  impurity  gradient  model  explains  the  formation  of 
p-type  regions  commonly  observed  at  the  interface  between 
an  n-type  epitaxial  layer  and  an  n^  substrate.  The  for¬ 
mation  of  heavier-doped  n-type  regions  located  at  the 
interface  between  an  n-type  epitaxial  layer  and  a 
chromium-doped  substrate  can  also  be  explained  by  the 
model . 

The  observed  screening  of  impurities  at  concentra¬ 
tions  above  the  intrinsic  carrier  concentration  is  also 
explained.  At  the  surface  the  surface  states  lower  the 
electron  concentration  and,  therefore,  increase  the  hole 
concentration.  This  increased  hole  concentration  screens 
the  impurities  at  the  surface  causing  the  experimentally 
observed  compensation  seen  in  Figure  4.3.  The  lower 
electron  concentration  at  the  surface  produces  band 
bending  and  an  electric  field  near  the  surface. 

This  electric  field  explains  the  chromium  pile-up 
observed  in  chromium-doped  epitaxial  layers  and  substrates. 
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The  model  predicts  that  this  electric  field  will  produce 
donor  and  acceptor  impurity  gradients  near  the  svirface. 
These  impurity  gradients  are  a  build-up  of  donors  and  a 
depletion  of  acceptors  at  the  surface.  The  electric 
field  at  the  surface  may  also  explain  why  GeJVs  epitaxial 
layers  are  self-compensated  by  acceptors,  which  are 
thought  to  be  donor  atom-vacancy  complexes.  Since  there 
is  a  depleted  concentration  of  acceptors  and  an  excess 
concentration  of  donors  at  the  surface,  from  mass  action 
laws  the  formation  of  donor-vacancy  complex  acceptors  will 
be  enhanced.  The  electric  field  will  cause  these  accep¬ 
tors  to  drift  away  from  the  surface  sustaining  the  deple¬ 
tion  needed  to  cause  enhanced  complex- formation. 

High  temperature  measurements  were  made  on  epitaxial 
layers  and  substrates  to  obtain  quantitative  results 
for  the  impurity  gradient  model.  A  self-consistent 
four-band  model  of  carrier  transport  was  used  to  calcu¬ 
late  an  intrinsic  carrier  concentration  of  5.8  x  lO^^cm  ^ 
at  1000  K.  Prom  this  four-band  model  heavily-doped  sub¬ 
strates  were  calculated  to  still  be  n'*'  and  chromium-doped 
substrates  were  calculated  to  be  slightly  p-type  at  1000  K. 

Using  these  data  p-type  interfacial  regions  500  to 

0 

1000  A  thick  are  predicted  for  the  growth  of  lightly- 
doped  n-type  epitaxial  layers  on  heavily-doped  substrates. 
Interfacial  regions  of  about  this  size  are  commonly 
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observed  at  the  n-n'*’  interfaces,  as  seen  in  Figure  1.1. 
This  impurity  gradient  model  should  be  applicable  to 
the  epitaxial  growth  of  materials  other  than  GaAs.  It 
should  be  noted  that  similar  0.5  ym  interfacial  regions 
have  been  observed  at  the  n-n^  interface  in  InP  [100] . 

By  combining  the  experimental  data  for  zero-bias 
capacitance  measurements  and  the  net  donor  profile  pre¬ 
dicted  by  the  impurity  gradient  model,  a  surface  state 
density  at  1000  K  was  calculated.  The  average  surface 

state  density  required  to  form  the  measured  impurity 
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gradients  was  5.1  x  10  cm 

This  calculated  surface  state  density  compares  fa¬ 
vorably  with  the  surface  state  density  of  6.3  x  lO^^cm” 
needed  to  explain  the  experimental  self-compensation 
data  of  Figure  4.3.  This  surprisingly  good  agreement 
between  two  different  experiments  tends  to  confirm  the 
impurity  gradient  model. 
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